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For  odd  degrees  of  freedom  the  characteristic 
function  of  a  Student-t  random  variable  is  expressible 
in  closed  form.   The  characteristic  function  of  an 
arbitrary  linear  combination  of  independent  t-variables 
is  then  derived  and  the  distribution  function  is  ob- 
tained, itself  expressible  as  a  weighted  sum  of 
Student-t  distribution  functions.   An  easy  method  of 
obtaining  the  weights  is  demonstrated. 

If  U-j^jU^,  .  .  .  ,Uj^  are  independent  random  variables 
and  X^  =  d+U^,  1  =  1,2,..., n  are  observable  random 


variables,  we  investigate  the  choice  of  a-,,ap,...,a   t 


o 


maximize  the  power  of  tests  of  the  form  a,  X-,  +  a^Xp  +  .  .  . 
+  a^X^  for  testing  Hq  :  d  =  0  against  H-,  :  d>0.   Some 
general  results  and  examples  are  given.   Of  particular 


vl 


Interest  Is  the  case  when  X.  Is  a  t-random  variable. 
One  application  is  in  a  two-stage  sampling  procedure 
to  solve  the  Behrens-Flsher  problem.   The  test  statistic 
has  the  distribution  of  a  weighted  sum  of  t-random 
variables.   It  is  shown  how  to  choose  the  weights  for 
maximum  power. 


CHAPTER  I 
INTRODUCTION 

1 . 1  Motivation:   The  Behrens-Flsher  Problem 

Given  that  the  distribution  of  the  population  JI. 

■  1 

2 
IS  normal  with  mean  y .  and  variance  a.  for  1=1,2,  let 

it  be  required  to  test  the  equality  of  y.,  and  \i^    based 

on  data  independently  drawn  from  n   and  n„.   If  we  know 

2 
the  value  of  6  =  °l/o^,    the  test  criterion  takes  the 


form 


(x-i  -x„) 
t  = ^— 


0   n,   np 


2    (n^-l)s^  +  eCn^-Ds^  _       p 

where  Sg  = ^  ^  ^  _  ^ and  n^,  x^  and  s^  are 

the  sample  size,  mean  and  variance,  respectively,  of 

the  sample  taken  from  H..   If  it  is  true  that  y.,  =  y^ 

then  t  will  be  distributed  according  to  the  Student-t 

probability  law  based  on  n,+np-2  degrees  of  freedom. 

Under  standard  analysis  of  variance  assumptions,  6  is 

taken  to  be  unity,  which  may  or  may  not  be  a  valid 

assumption.   If  not,  the  resulting  test  will  be  biased 

to  a  degree  commensurable  with  the  deviation  from  that 

assumption.   When  the  value  of  6  is  unknown,  the  test 

criterion  has  an  unknown  distribution.   This  situation 


is  commonly  referred  to  as  the  Behrens-Fisher  problem. 
Under  Behrens-Fisher  conditions  no  criterion  with 
distribution  as  simple  as  the  t  above  is  available. 
Fisher  proposed  a  fiducial  test  based  on  the  so-called 
Behrens-Fisher  statistic, 

(x   -  Xp) 

°  ^  -^ ^-U2  (1.1.1) 

!l  +  !2 
n^  ""  "2 

and  Sukhatme  [16]  calculated  approximate  percentage 
points.   Other  solutions  have  been  proposed,  notably  by 
Welch  [17],  Scheffe  [13],  Welch-Aspin  [l8],  McCullough, 
et  al.  [9],  and  Cochran  [3],  to  name  a  few.   Many  of 
the  proposed  solutions  are  based  on  an  approximation  or 
asymptotic  expansion  of  the  distribution  of  the  statistic 
D,   Fisher  rev;rote  this  quantity  as 


(1.1.2) 


random  variables.   Ray  and  Pitman  [11]  developed  an 
explicit  form  for  the  density  function  of  D.   Using  a 
series  expansion  for  the  density  of  a  weighted  sum  of 
chi-square  random  variables,  they  produced  a  very 
complex  formula  representing  the  density  of  D  and 

showed  how  one  might  compute  percentage  points.   Using 

2       ? 
the  fact  that,  conditional  on  s   and  s^,  D  is  the 


weighted  difference  of  two  Independent  t-varlables, 
Rahman  and  Saleh  [10]  expressed  the  density  In  terms 
of  Intergrals  of  Appell  functions  and  used  numerical 
quadrature  techniques  to  table  critical  values.   In 
each  of  these  attempts  an  approximation  was  used  some- 
where In  the  development  before  any  values  could  be 
tabulated.   In  addition,  It  is  a  very  complex  procedure 
to  obtain  non-tabulated  critical  values  using  any  of 
these  methods. 

Although  It  is  still  not  known  whether  there  exists 
a  fixed  sample  size  test  whose  significance  level  is 
Independent  of  the  unknown  parameters,  Dantzlg  [i|] 
proved  there  exists  no  fixed  sample  size  test  whose 
power  is  independent  of  a^.   A  two-stage  sampling 
scheme  vias   devised  by  Stein  [151,  the  power  of  the 
resulting  t-test  being  independent  of  the  nuisance 
parameters.   This  approach  may  be  modified  to  give  a 
solution  to  the  Behrens-Plsher  problem  (see  Sections  1.2 
and  6.1),  and  in  so  doing,  the  test  criterion  has  the 
distribution  of  a  weighted  difference  of  Student-t  vari- 
ables.  Chapman  [2]  tabulated  some  values  of  the  cumula- 
tive distribution  function  for  this  weighted  difference 
but  only  for  the  special  case  of  equal  weights  and  equal 
degrees  of  freedom.   Ghosh  [71  also  investigated  this 
distribution  for  the  case  of  equal  weights.   He  gave 
explicit  formulas  for  finding  the  distribution  function 


when  the  degrees  of  freedom  are  equal  and  less  than  5, 
and  tabulated  a  few  values  of  the  probability  integral 
when  the  degrees  of  freedom  are  equal  and  greater  than 
5.   In  addition  he  used  a  series  expansion  to  obtain  an 
approximation  to  the  cumulative  distribution  function 
when  the  degrees  of.  freedom,  differ.   Ghosh,  incidentally,, 
pointed  out  that  many  of  Chapman's  values  were  incorrect. 

Thus  far,  the  distribution  of  a  difference  of  two 
t-variables  with  different  weights  has  not  been  investi- 
gated in  the  literature.   Furthermore,  if  the  degrees  of 
freedom  associated  v;ith  the  two  t-variables  differ,  the 
exact  distribution  has  not  been  derived.   The  Behrens- 
Fisher  problem  would  be  solvable  (section  1.2)  if  the 
distribution  of  a  v/eighted  difference  of  Independent 
t-varlables  with  arbitrary  weights  and  degrees  of 
freedom  vjere  known. 

1 . 2   A  Sequential  Solution  to  the  Behrens-Fisher  Problem 

Given  two  normal  populations,  IT,  and  !!„,  where  the 
mean  and  variance  of  IT.  are  y.  and  of   respectively  (1  =  1,2), 
we  desire  to  test  the  hypothesis 

To  this  end  one  might   take    n.  observations  from  IT. 
and  use  the  statistic  D  In  (1.1.1)  as  test  criterion. 
If  the  sample  sizes  are  small  however,  the  usual  normal 
approximation  may  not  be  accurate. 

Consider  the  following  two-stage  sampling  scheme: 


sample  varlates  (X^.},  j=l, 2, . . . , n^ ,  from  n^ ,  1=1,2, 
and  compute 


h  -  I 


X.  . 
J=l   n^ 

n. 


T-r        n2 


ss.  =  y  (X.  .-X.  )■ 


J.  =  max  n.  +  1  , 
1        1    ' 


SS. 
1 

2  2 
a.  T 
_  1 


+  1 


v;here  a^,    a^  and  x  are  preassigned  constants  and  the 
square  brackets  indicate  "the  integer  part  of." 

The  second  stage  consists  of  sampling  the  additional 
observations  {X^.},  j=n^+  1  , n^ +  2  , . . . , M^  f rom  n .  (i  =  l,2). 
The  test  criterion  for  testing  H„  will  be 


T  = 


a^ .X^  . 

2j  2j 


(1.2.1) 


where   'tct^.},  j=l,  2  ,  ,  .  .  ,N^  ,  are  constants  satisfying 


a.  -,  =a. 


il   12  • 

J. 

1 

ii).   I_a,  .  =  1 

N, 


N. 

1 


J=l  ^J 


(1.2.2) 


ill).  I    a^ 


j=l 


SS. 

where  f.  =   -^=; — ^ 
1      2  2 
1  a.  T 


Such  constants  may  always  be  found  since  the  minimal 

^^i        1  11 

value  of   V   2   is   i   and  N.>f.  so  that  ^<^   . 
_I_^a..      N.        1   1         N.  f . 

Letting  d=  \i-.-v^p   it  is  shown  in  section  6.1  that  T 


may  be  represented  as 


where  ^j^=  J^j^- 1  and  t.  (1  =  1,2)  are  Independent  t  random 
variables  on  v^  degrees  of  freedom.   Thus,  when  Hq :  d=0 
Is  true,  the  test  criterion  has  the  distribution  of  a 
weighted  sum  of  t-varlables.   Letting  X  (a.,  ,  ap,  n,  ,  n^  ) 
be  the  100(l-(j))  percentage  point  of  the  distribution  of 
T,  we  reject  H^  In  favor  of  the  one-sided  alternative, 
H^  :  d>.0,  whenever  T>X  ,  (a.^  ,ap  ,n-,  jn^  ) .   The  experimenter 
Is  free  to  select  the  constants  {a..},  j=l, 2 , . . . ,N. ,  a. 
and  n.  for  1=1,2,  as  well  as  x.   The  Initial  sample 
sizes,  n. ,  Influence  the  expected  total  sample  size, 
^(N-j^+Np)  ,  while  T  and  the  a .  '  s  Influence  the  power 
function  of  the  test  T.   A  discussion  of  a  discreet 
choice  of  T,  n,,  n^,  and  the  a..'s  Is  deferred  until 
section  6.1.   Presently  we  are  concerned  with  the  selec- 
tion of  a-,  and  ap.   Two  problems  of  Immediate  concern, 
then,  are 

1.   that  of  obtaining  the  percentage  points,  A,, 
of  the  distribution  of  T,  and 

11.   the  method  of  selecting  a.and  ap  to  maximize 
the  power  of  T. 

1.3.   Statement  of  the  Problem 

Let  w, ,Wp,...,w  be  an  arbitrary  set  of  constants 
and,  for  1=1, 2,..., n,  let  t.  be  Student-t  random  variables 
Independently  distributed  on  v.  degrees  of  freedom.   For 


d>0,  denote  by  T   the  linear  combination 

^d  "  w^(t;j_+d)  +  W2(t2+d)  +.  .  .+  w^(t^+d)     (1.3.1) 

which  may  be  written  as 

n 
T^  =  U  +  d  ^  w.  ,       •   :        (1.3.2) 
1  =  1 
n 
where  U  =  I   vi.t..      The  first  problem  above  may  be 
1  =  1  ^  ^ 

generalized  to  that  of  finding  the  distribution  of  T , . 

The  second  problem  becomes  that  of  selecting  weights, 

{w.  }j  1=1, 2,..., nj  to  maximize  the  power  of  the  test  T, 

for  testing 

Hq:  d=0 

vs.    H-^:  d>0  . 

The  purpose  of  this  current  work  is  to  investigate 

the  solutions  to  these  problems.   VJe  show  in  particular 

that  the  probability  density  of  T,  may  be  written  as  a 

linear  combination  of  Student-t  density  functions  whose 

weights  depend  on  {w. } , 1=1, 2, . . . ,n,  { v . } , 1=1, 2 , . . . , n, 

and  n.   The  power  of  T,  depends  also  on  these  quantities. 

If  v.=Vp=...=v   it  is  shown  that  no  linear  combination 

of  the  t.'s  has  uniformly  greater  poiver  than  for  the 

case  when  Wt=w„=...=w  .   If  v.j^v.    for  at  least  one  pair, 
12      n       1  .  J  t^        y 

i?^j  ,  no  distinct  set  of  weights,  {w.},  exists  which  will 

maximize  the  power  of  T,  for  all  size  tests. 
^  d 

Chapter  2  is  an  investigation  of  the  properties  of 

the  distribution  of  t .  .    ,  including  derivation  of  the 

1 
characteristic  function  and  distribution  function  in 


closed  form  for  odd  v..   Chapter  3  includes  the  deriva-- 

tion  of  the  distribution  of  T,  for  odd  v.;  the  case  of 

d  1 

even  degrees  of  freedom  is  considered  separately.   For 
the  case  n=2,  tables  of  the  percentage  points  of  T,  are 
given.   Some  general  theorems  for  comparing  the  power  of 
two  tests  are  given  in  Chapter  ^  along  with  some  applica- 
tions, including  several  examples  and  counterexamples. 

In  particular  it  is  shown  that  t.,    has  greater  Dov/er 

1 

than  t.,   v.'henever  v .  >v  .  .   Chapter  5  gives  a  comparison 
J<^j  1   J       ^       ^ 

of  results  with  the  normal  approximation  and  with  Ghosh's 
results.   Finally,  along  with  the  solution  to  the 
Behrens-Fisher  problem,  another  application  is  indicated 
in  Chapter  6. 

1.4   Notation 

The  following  notation  will  be  used  whenever  con- 
venient : 

a.  Abbreviations 

cdf  cumulative  distribution  function 

d..  f.  degrees  of  freedom 

lid  independent,  identically  distributed 

iff  if  and  only  if 

wp  p  with  probability  p 

r.v.  random  variable 

b.  Symbols 
(A). .  the  element  of  row  i,  column  j  of 


ij 


matrix  A 


(A)  the  i    row  of  matrix  A 

th 

(A)  .  the  j    column  of  matrix  A 

•J 

EiX)  expected  value  of  X 

7(X)  variance  of  X 

Pr{A}  the  probability  of  event  A 

r{x}  the  gamma  function 

B(x,y)  the  beta  function,  ^[{^+  ) 

J{X^Y}  the  jacobian  of  the  transformation 
from  X  to  Y 

*         -1 

F  F   ,  the  Inverse  function  of  F 

C^(P)  the  coefficient  of  x^  in  P 

[k]  reference  k  in  bibliography 

[k:l]  page  or  line  1  in  reference  k 

=^  implies 

<=>  is  equivalent  to 

~  is  distributed  as      .  ■ 

V  for  all 

5  such  that 

therefore 

§  section 

Generic  random  variables  and  their  distributions 

t  Student-t  r.v.  on  v  d.f. 

V 

X^  pre-t  r.v.  (  v/^  ) 

2 

X  chi-square  r.v.  on  v  d.f. 

Z  standard  normal  r.v. 


density  of  Xp  .  , -, 


h. 
J 

H.         cdf  of  Xo  .^. 
J  2j+l 


CHAPTER  2 
THE  PRE-t  DISTRIBUTION 

2. 1   Density  and  Moments 

We  denote  by  t^  a  random  variable  (r.v.)  having 
the  Student-t  distribution  on  v  degrees  of  freedom 


(d.f.).   The  density  of  t   is 

V 


f(t)  =    y  '-  (1+  ^ 

Making  the  transformation 

t 

y^   =  —  , 


v+l 
2     2 


we  obtain 


Jit^-^X^}   =    AJ  , 


so  the  density  of  X   Is 

V 


v  +  l 

2 


2'2 

The  r.v.  X^  shall  be  referred  to  as  the  "pre-t  r.v. 
on  V  d.f."   We  use  this  in  place  of  the  t  In  the 
ensuing  theoretical  development  to  simplify  the  math- 
ematical formulations.   It  is  clear  that  h(x)  is 
symmetric  about  x=0  and  decreasing  as  |x|  increases. 
Hence,  all  odd  moments  of  X   are  zero.   For  k=l,2,3,,. 

10 


11 


2'   '2 
v>2k         (2.1.1) 


In  particular. 


^(X  )  =  0 

V 

and  V(X    )  =  -^  ,   v>2. 
V    v-2  ' 

2 
If  Z  represents  a  standard  normal  r.v.  and  x  > 

an  Independent  chi-square  r.v.  with  v  d.f.,  we  may 

Z  2        2 

represent  X^^  by  —— -.   If  v„  >  v,  and  y    and  y 

"J  X^  2    1      '^v-^      ^2~^1 

denote  chi-square  r.v.'s  independently  distributed  on 
v^  and  '^2~^l  '^■^•^    respectively,  then 

■7       2.2      2      ~^ 

^  >  Z  (x    +  X,,  _,,  )    with  probability  1. 

x;.       ""l    ''2  ''l 

1 

Thus  for  any  t>0,   ' 

Pr{|X^  |>t}  >  Pr{|X^  |>t}  whenever  v  <v  .     (2.1.2) 

That  is,  Ix   lis  stochastically  less  than  Ix   |. 

The  remainder  of  this  chapter  will  be  devoted  to 

deriving  the  characteristic  function  and  cdf  of  X   and 

exploring  other  useful  properties  of  its  distribution. 

2.2   The  Characteristic  Function  of  X 

. _ . ^J 

We  may  use  a  property  of  the  Bessel  functions  to 
give  a  closed  form  representation  of  the  characteristic 
function  of  X   in  the  case  when  v  is  odd.   The  modified 

V 

sperical  Eessel  function  of  the  third  kind  is  given 
from  [l:9.6.25]as 


12 


K^(xz) 


r{f+|}(2z) 


cos(xt) 

.,2  ,  ,2.f+ig  ^^ 
0  Ct  +  z  J 


for  f  2  -2  and  |arg  z  |  <  |-  .   Letting  f  =  ^   (v  >  -l) 


and  z  =  l,  x= I  6  I ,  we  obtain 


v+1.  .^/2 


K^  (|e|)  = 


r{^}  2 


so  that 


TT  -  e 


v/2 


cos(  |e|t)   ,, 
^0   (l+t2)^-^(--^l)  '"' 


(2.2.1) 


CO 

0   ^^2"'2"'' 


%I<^V'^" 


(2.2.2) 


The  characteristic  function  of  X   is 

V 

(})(e;v)  =  ^(e^^^^) 

=  ff(coseX  +  i  sinex  ) 

V  V 

=  ^(cosex  ), 


since  ^(sinex 


sln0x'h(x)dx  =  0,  h(-)  beln^ 


symmetric  about  zero  and  sin(0x)  an  odd  function.   We 
have,  therefore. 


(e;v)  = 


cosex  .^^  2  -'"^^^+1) 

(1+x  )       dx 


-co  3(2'2-' 
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rOO 


=    2 


cos(xl  e  I  )    ,^^    2 


-^-sCv+l) 


0      6(i,| 


(1+x    )  dx    , 


or,    from    (2.2.2), 


(e;v)  =  ^(^4^1  "  K  del)   . 


(2.2.3: 


By    [1:    10.2.15,    10.2.11] 


-fel      V;,^^)    =   -TF      ^"'    I    (n+^^„k)(2z)-^ 


n 
k=0 


(n+k) ! 


where    {n+h,k)    =   y)  (y^_{A  i     y    k=0  , 1,  2,  .  .  .  ,  n.       Substituting, 


r-  I 


(n+k)!       ,o„^-k 


V..^(z)    -   ^      e--^£Q    k!(n-k)! 


(2z)    "    , 


or, 


v..(^)  =1^1  -'\l^ 


z    V       (n+k) !       ^o„^-k 
^k! (n-k) ! 


(2z)    '\    n=0,l,2,.. . 
(2.2.4) 


Letting  v = 2ra+l  for  m=0,l,2,...,  and  using  the  identity 


{2m±l}  2'"'^^'"^^'*  =    (2m)!/2¥ 


m!  2' 


(2.2.3)  becomes 


(0;2m+l)  = 


i M 


,^2m+ljl  2 


J'2(2m+1) 


•^m+%<l«l> 
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ig(2m+l) 


2  e 


k=0 


(m+k) ! 
k! (m-k) ! 


(2|0|) 


-k 


(2|9|)"^  m!   -I  6  I  ■?   (m+k)!   .pi.^-k 
(2m)!    ^     ^£q  k!(m-k)!^2|^l^ 


e-|e|  m!    '   (y)!   (2|e|)' 
(2m)  !i^^Q  k  !  (m-k)  !  ^  '  '  ^ 


or  finally. 


(e;2m+l)  =  e' 


IliOTj^^Q  k!(m-k)! 


(2|e|)^^  .    (2.2.5) 


m 


Denoting   q  (z) 
simply  write 


m!         V    _(2m-k)!     ^^^^k 
k  =  0 


m-'         r2iyr,  L-l^TTiPkyi^^z)    ,    we   may 


(e;2m+l)    =    e    '     'q^( |e|)    ,    m=0,l,2, 


In  particular,    v/e   have 


(2.2.6) 


;1)  =    e" 

i;3)  =    e' 

;5)  =    e' 

;7)  =   e" 

;9)  =   e" 


(1  + 
(1  + 
(1  + 
(1  + 


el) 


4le|^) 


+  f|e|^+I5le|3) 


■^  I      I  2  2 


3+  i_ 
105 


(2.2.7) 
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2. 3   The  Cumulative  Distribution  Function 
After  an  Integration  by  parts. 


cos^^'e  de  =  4  cos^^^-^esine  +  ^"^-^ 


2m 


2m 


2m-2n  ,. 
cos     9  de 


and. 
Therefore, 


2        1 
cos  6  de  =  2^sln8cose  + 


cos2%    de    =   1^  cos2"^-leslne    +    2m{Tm^-2)    ^os^^^'^eslne 


,         (2m-l)(2m-3)    ^^    2m-5n    •    n    ^ 
^    2m(2m-2)(2m-H)    '^^^  ^^^^^^    +    • 


+      (2m-l)(2m-3)- • •(5)(3) 
2m(2m-2)(2m-^)-  •  •  (  6  )  ( ij ) 


2 
cos    e    d( 


=    sine.         2m-l  2m-l  2m-3n     ,     (2m-l)(2m-3)  2m-5, 

2m    ^^°^  ^    +    2^32    ^°^  0    +    (2m-2)(2m-4)    °°^  * 

+  +    (2m-l)(2m-3)---(5)(3)      ,^    „. 

•••    ^(2m-2)(2m-l|)...(6)(M(2)°°^^^    . 


+    (2m-l)(2m-3)--'(5)(3) 
2m(2m-2). . . (4)(2) 


(2.3.1) 


By  Gauss  duplication. 


r{m+i,}  =  il5ll2m} 


2  2'^r{m} 


(2.3.2) 


so  that 


,1   2m+l^ 
'^2'~2~ 


r{m+l} 


r{i}r{^5^}       /^  rfm+Jg} 
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2m  {(m-l)!2'^  ^} 
i\  (2m-l)  ! 


2rn 


(2m-2)(2m-4)---(6)(H)(2) 
(2m-l)(2m-3)---(5)(3)(l) 


(2.3.3) 


Denote  the  density  of  X„  ,,  by  h   and  Its  cdf  by  H  , 

2m+l   '^   m  m 


m=0,l,2. 


For  t  >  0  and  9  =  tan   t 


H^(t)  =  ^  + 


-1^ 


1 


1  2m+l 


p  -(m+1) 
(1  +x  )       dx 


0  ^'2'      2 

'  1 
ri  2m+l. 
4'  2  ^ 


2in 
cos   u  du 


i+  i 


2m(2m-2)(.2m-J^)-  •  •  (  6  )  ( i) )  (  2  ) 


(2m-l)(2m-3)---(5)(3)(l) 


X  I 


where  I  represents  the  right-hand  side  of  equation  (2.3.1) 

The  closed  form  expression  for  the  cdf  of  X^  ,,  Is 

2m+l 

therefore 


H^(t)  -  ^  +  - 
m      2    TT 


+  sin  6  {cos  6  +  Icos^e  +  i^li^P-cos^Q 
3        (3) (5) 


+     +  (2)(4)---  (2m-4)(2m-2)„_2m-l„, 
^  •■•  ^  (3)(5)---(2m-3)(2m-l)^°^     ^> 


where,   9  =  tan  t,   t>0. 


(2.3.M 


2 . 4   Additional  Properties 

The  following  properties  about  the  distribution  of 


X  will  be  useful  In  later  sections 


Property  1. 


9h(t) 


t(v+l) 
3(2->2^ 


(1+t^) 


■Js(v+3) 


Property  2.    H  ,T(t)-H  (t)  =  75-^  h  (t),  m=0,l,2, 
m+1     m      2m+l   m    '     35. 

t>0 
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Property  3 


Hjt)-Hj^(t) 


m-l  h. (t) 


=  t  I 


2J+1 


m  >  .k. 


Property  ^ 


(e;2m+l 
39  . 


2m-l  ,    '  ni-l>2,3. 


2 
Property  5.   .cj)(0;2m+3)  =  0(e;2m+l)  +  (2m-l^J2m+lj 

m=l,2,3, . . .  - 
Property  1  follows  by  differentiating  both  sides  of 

h(t)  =  —Y^    (1  +  t  ) 
3(2"5  2''' 

To  prove  Property  2  xve  may  use  the  results  of 
section  2.3: 


H  ^T (t)-H  (t) 
m+1     m 


_  1_.^J  (2)(4)(6)--  -(2111)  „„  2m+l 


(l)(3)(5)---(2m+l) 


cos     9 


(9  =  tan  H) 


T  ,^,^111x2        -(m+1) 
TT  (2m+l )  ! 


^  r{m+l}    2^"^r{m}  (i  +  t^)"^'""'^^ 


^"""^^  r{|}   r{2m}2v/? 


Using  (2.3.2),  we  obtain 


H  ^. (t)-H  (t 
m+1'     m 


t      r(m-H)     n^.f^l"'"'"'^' 


TT^  h    (t)  ,  m=0,l,2,  .  .  . 
2m+l   m     '     }    y    } 


If  m  >  k. 


m-1 


H  (t)-H  (t)-  =   J  {H  ^,   .(t)-H  ^,   .  ,(t)} 


j-k 


m-l 


J,   ^2m+2k-2.1-l  \+k-,1-l^^^^ 


j=k 
m-l 


1 


which  Is  Property  3.   To  demonstrate  Property  ^,  we 
may  begin  by  writing,  from  (2.2.6), 

^hL^D  o   n^    (m-l)!  "^v^  (2m-k-2)!  /o,.|^k 
^   'Me;2m-1)  =    ^g^_g^,  J^  ^;(^_i,_i);(2|9|)   , 

a  polynomial  of  degree  m-l  in  |9|.   Also, 


/Q-i_,T>     -|9|  in!     V    (2m-k)!  ,„|.|vk 
^^'2^+1^  =  ^     T2i^  ,^,  k  (m-k)!  ^2|e|)  , 

k=0 


so  that 


9|  3c^(e;2mH-l)  ^   m!    ^  ^2m^  ^^k  |  9  |  ^"1- I  9  |  ^  , 
3  9         (2ra)!  ,^„  k!(m-k)!    ^  i  i     i  i  ^> 
k=0 


a  polynomial  of  degree  m  in  | 9 | .   Using  C      to  denote 

-p 
"the  coefficient  of  x   in",  we  have       •  - 


^r   ^9(j)(e;2m+l)^  ^  g- 


9  I  m! 
(2m)! 


(2m-r-l)!      r+1.    ^ 
+l)!(m-r-l)!  ^    ^^^^^ 


(2m-r) !2- 
r ! (m-r ) ! 
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m!    (2m-r-l)!2^ 


'    '  '(-2^  t^q^iyrtf^f:^(2(r4.1)(m-r)-(2m-r)(r+l)} 


-lei  m! (2m-r-l) !2^ 


=  -^     (2m)!(r-l)!(m-r)!  '  r=l , 2 , 3 , . . . ,m  .   (2.4.1) 
Also, 


^^   r      9   w  ^ .  -^._  -.v-,        1    -lei  (m-1)! 
^lel^-  -(2¥rn*^^''^'-^^^  =  -  T2i^Ty^  '  '(2m-2)! 


(2m-r-l)!   2^-^ 


(r-1) ! (m-r) ! 


-  ^-|9|   m!(2m-r-l)!  2^      _^  „ 

-  -^     (2m) ! (r-1) ! (m-r)!  '   ^"1.2, ...,m 

which  is  identical  to  (2.4.1).   When  r=0,  we  have 


0  r    9(t.(e;2m+l)i  ^   -|e|   m!   (2m-l)!2    (2m)! 
leT     90  (2m)!   (m-1)!   "   (m)  ! 


=  0 


which  proves  Property  4. 

Property  5  may  be  proved  by  induction.   From 
(2.2.7)  we  have 


(e;5)  =   e"l^l(i  +  |e|  +||e|2) 


(e;3)  =  e-i"i (1  +  |e| ) 
(e;l)  =  e-|0|   ,     ■■ 
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so  that . 

<|)(9;5)  =  4>(e;3)  +  -^yjj^ycj)  ( 6  ;  1 )  , 

which  shows  that  Property  5  holds  for  m=l.   The 
induction  hypothesis  is 

c|)(e;2m+l)  =  (f,(e;2m-l)  +  (  2^-3  )  (  2m-l)'^  ^  ^  ^  ^"^"^  ^  ' 

(2.  A. 2) 

Writing 

0(9;2m-l)  =  ||^±lj(|,(e;2m-l)  -  ^^-|_^  (  9  .  2m-l )  , 
(2.4.2)  becomes  ' 


<|)(9;2m+l)  _    1   ,  ( p,  .  Or.    i\    +       6^       ^(ej2m-3) 
(2m+l)      T2Hriy*^9'2m-l)  +  (2^_i)(2^+i)    (2m-3) 


~    -  T2irTyr2iTry*^e'2m-i)  .       (2.1J.3) 

From  Property  4,     . 


fl2 
(J)0;2m+3)  =-4.(e;2m+l)  +  (2m-l)  (2m+l)'^^^ '^''-^' 

if  and  only  if  for  all  8^ 


9(i)(9;2m+l)  ^  9  4)(9;2m-l)         9^       (t)(9;2m-3) 
2m+l         2m-l      '  (2m-l)(2m+l)   (2m-3) 


T^ -I  N  /  ^ — rT~v<)^  (  9  j  2m-l ) 

(2m-l) ( 2m+l)    ' 
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which  is  equivalent  to  (2.i|.3).   Thus  Property  5  is 
established. 

Property  H    is  useful  in  demonstrating  the  infinite 
divisibility  of  the  distribution  of  X  .   A  function  V 
is  the  characteristic  function  of  an  infinitely 
divisible  distribution  if  and  only  if  'l'(6)=e~"    , 
where  to  has  completely  monotone  derivative  and  a)(0)  =  0 
(see  [6:  p.  4251).   Grosswald  [8]  has  shown  that 

64'(e;2m-l)    '   -n  i  ?  '   ' 

ct.(e;2m+l)   '   ni-0,1,^,... 

is  a  completely  monotone  function  of  9.   Using  Property 
^,  we  must  have 

9(i)(6;2m+l) 
(|)(e-2m+l) "  g^(-log  (j)(e;2m+l)  )   is  completely 

monotone  in  9.   That  is,  -log  (j)  (9  ;  2m+l )  has  completely 
monotone  derivative.   Thus, 


^(9)  =  e^°S*^^'2m+l)^  cD(9;2m+l) 


is  the  characteristic  function  of  an  infinitely  divisible 
distribution.   That  is,  the  pre-t  r.v.  with  odd  degrees 
of  freedom  has  an  infinitely  divisible  distribution. 


CHAPTER  3 

THE  DISTRIBUTION  OF 
A  LINEAR  COMBINATION  OP  t-VARIABLES 


3 • 1   Simplification 

Letting  t^  denote  a  t  r.v.  on  v  d.f.,  the  Cauchy 
variable  may  be  represented  by  t,  and  the  standard 
normal  by  t^.   For  an  arbitrary  set  of  constants, 
c^,C2,...c^,  the  linear  combination 

n 
T,  =   y  c.t,  . 

has  the  Cauchy  distribution,  and. 


n 

T^ 

1  =  1 


T  =    y  c .  t 

2    .  '--,  1  ^ 


has   the   normal   distribution,    vihere   t    .    are    lid   as    t    , 

VI  v' 

1  =  1, 2,..., n.   If  1  <  V  <  0°  ,  no  such  simple  distribution 

n 
exists  for  the  linear  combination   yet.. 

1=1  ^  ^^ 
Here  we  consider  the  generalized  linear  combination 


n 
T'  =   y  c.t 


where  the  t's  are  Independently,  although  not  necessar- 
ily Identically  distributed.   The  cumulative  distribu- 
tion function  of  T'  is       . 
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n  n 

Pr{  I   c    t       <yn    =   ?v[    I   c.A^.X^^<y'} 
1=1    1  1=1  1   1  ^1 


(3.1.1) 


where  X   denote  independent  pre-t  r.v.'s  on  v.  d.f, 

1 


Letting  oj  =   T  C-A;".  ,  a. 

'■'Ill 

1  =  1     ^     ^ 


c  .  /v . 
1   1 


and  y 


LL 


(3.1.1)  becomes   Pr{  I    a^X^  <  y}.   For  m.  =0,1,2,... 
1  =  1     1  "'" 

and  V.  =2m.+l,  the  characteristic  function  of  a.X 

11  IV. 

1 
Is,  from  (2.2.6), 

-  I  a .  e  I 
(^(a.e;v.)  =  e      q^(|a.e|)  , 

so  that  ct.X    and  |  a .  |  X    have  identical  distributions 

1  i 

We  therefore  only  need  to  consider  the  distribution  of 


T  =   y  a.X 


J-1  '    ^J 


(3.1.2) 


for  the  case  a.  >  0,  j=l, 2, . . . ,n,  and  I    a .  = 1. 
•^     .  j=l  -^ 


3.2   The  Distribution  of  T  Using  the  Q  Matrix 
The  characteristic  function  of  T  is 


leT 


(6)  =  Ee^^    =   E.   expde  I    a   X      } 

j=i  -i   ""j 


n    iBa.X^. 
=   n  ffe   J  ^J 

J=l 


n  d)(ea.  :  V. 
J=l     J  ' 


V/ith    V    =  2m.+l    (m  .  =  0, 1,  2  ,  .  .  .  )  ,     (2.2,6)    gives    us 


n     -\ea.\  -loin 

.(6)  -  n  e     J  q   (|6a  I)  =  e 


2ii 


n  q^  (I ea, I )  . 


Denote  by  $_  and  0(x)  the  arrays 
'^(e;l) 

4>(6;5) 


J-1      J 


1 
Gx 

ex 


m,.     '       J 


(3.2.1) 


and   0   -    0(1).       Let    $   be   a   matrix   v;hose    element    in   row   i, 

,   -|e| 

column  j  is  the  coefficient  of  |e|^e      in  (!)(G;2i  +  l) 

(or  simply  the  coefficient  of  j  6  j  •^"  in  q.(|e|)  ),  for 
j=0,l,2,...;  1=0,1,2,....   We  may  now  rewrite  (2.2.6)  as 


QO 


(3.2.2 


th 


Denoting  the  1*^^'  row  of  Q   by  (Q)  .     ,  we  have 
Qpdel)  =  ^^)r.^  and  q   (|ea  I)  =  (Q)         0(a), 
so  that  (3.2.1)  is 


-|e|  n 
*T^Q^  =  e    '  n  (Q)        0(a)  .    (3.2.3) 

j=i       y      '^ 

|9| 
Clearly,  e   <Prp(Q)    Is  a  polynomial  in  |e|,  so  we  may 

find  a  vector  i   such  that 


(3.2.li) 


The  matrix  Q   is  lower  triangular  and  is  non-singular, 
Now,  (3.2.2)  becomes  Iqi 


Q    -^^    =    e 


0 
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so  that 


(p^iQ)    =    e     £'0  =  VQ    ^$ 
or  finally, 

$,p(0)  =  n'i  (3.2.5) 

where  n'  =    V  Q~    .   Thus,  we  may  v/rlte  the  characteristic 
function  of  a  linear  combination  of  independent  t's  as 
a  linear  combination  of  Student-t  characteristic  functions 
Since  (?)     has  zeroes  after  the  first  m.+l  elements 
(the  first  elem.ent,  corresponding  to  column  0,  is  1), 


(S)    0(a.)  is  a  polynomial  of  degree  at  most  m.  in  |e| 
j  •    J  J 

n 

^  ^^^m   ®(o'n-)>  therefore,  is  a  polynomial  in  \d\    of 
j=l   ™j'    -^ 

n 
degree  at  most  s  =   J  m.  .   This  tells  us  that  £  has  at 

j=l  J 
most  s  non-zero  elements  after  the  first  element  (which 

is  1).   Therefore,  n  is  a  vector  with  zeroes  after  its 

first  s+1  elements. 

The  cdf  of  T  may  now  be  written  as  the  linear 

combination  of  pre-t  cdf's: 

H(y)  =  Pr{T  <  y} 

=  ngPrfX-^  <  y}  +  Ti^Pr{X^  <  y}  +  .  .  .  +  ngPr{X2g^^  <  y} 

or  simply,  using  the  notation  of  section  2.3, 


H(y)  =   I  n.H„.,^(y)         (3.2.6) 
Setting  G  =  0  in  (3.2.5)  or  y  =  0  in  (3.2.6),  we 
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see  that   )^  n .  =  1 .   In  addition,  some  preliminary 

1  =  0 
results  and  a  computer  investigation  indicate  that 

n .  is  always  non-negative,  i  =  0 , 1 , 2  ,  . . . , s . 

The  first  eleven  rows  and  columns  of  Q   and 

Q~     are  given  in  Tables  3-1  and  3-2  for  convenience 

in  calculating  the  n^'s.   An  easy  way  to  obtain  the 

elements  of  Q        is  presented  in  the  next  section. 


3 .  3   The 


Matrix 


In  order  to  obtain  the  elements  of  Q~      one  may, 
of  course,  proceed  by  solving  QxQ~     =1,  the  identity 
matrix;  Q    being  lower-triangular,  this  is  easily 
accomplished  row  by  rov;  in  a  sequential  manner.   A 
more  general  method  is  given  here. 

For  m  and  j  in  {0,1,2,...  }  let  §  be  a  matrix 
whose  element  in  row  m  and  column  j  is  given  by 

(-l)^"J(2j)!(m+l)! 


^^'Vj   =  ' 


2'"(m-j)!J!(2j-m+l)! 


0  ,   otherwise 


,  J  <m  <  2j+l 


(3.3.1) 


*    -1 
It  will  be  shown  that  Q     =  Q      . 

The  (m-k)    element  of  Q   is  the  coefficient  of 

|e|*^  in  q^(l6|),  for  m,k  in  {0,1,2,...  ,j.   That  is. 


(Q) 


m,k 


m!(2m-k)!  2^ 
(2m) !k! (m-k) ! 

0  ,   otherwise 


m  >  k 


(3.3.2) 
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The  j    rov;  of  Q   is  given  by 

r 


^^h-  -xkr 


(2j)!n0   (2j-l)!-,l   (2J-2)!„2 


2' 


0!j!"  'l!(j-l)r  '2!(j-2)!^  '••'' 


th  * 

The  k    column  of  Q      is 


j!  oJ 


j  !0! 


2^  0,  0,  0, 


(3.3.3) 


(2k)! 


k     k! 


0 


(-l)°(k+l)! 
2^^0!  (k+1)  ! 

(-l)^(k+2)! 
2^+hlk! 

(-l)^(k+3)! 
2^^^2!(k-l)! 


(-l)^^^(2k+2)! 
2^^^'^^(k+l)!0! 


k  zeroes 


(3.3.M 


Clearly^ 


(e)j.x(5'),^ 


0  If  J  <  k 

1  if  j  =  k 
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(3.3.5) 


For  j  >  k  denote  ( Q) .  x ( $  )  ^  by  ^ .  If  k  <  J  <  2k+l, 


.  =  (2k)!j!'^\^.    i   (2j-k-l).!(k+l  +  l 
"-        k!(2j)!^^Q^  ^^   i!(j-k-l)!(k-l+l 


and  If  j  >  2k+l, 


f  =  (2k) !j!  ^y^.  ^.1   (2j-k-i)!(k+i+l) 
^    k!(2j)!  .^Q^~^^   l!(j-k-l)!(k-i+l) 


or  more  concisely. 


r  =  (2k)!j!   ?  ,_  .i   (2j-k-l)!(k+i+l) 

*-    k!(2j)!  .£q^  -^^   l!(j-k-l)!(k-i  +  l)!  '  (3.3.6) 


where  5  =  ■ 


j-k   if   k  <  j  <  2k+l 
k+1   if   j  >  2k+l 


Surely,  q'     =    Q   ^    iff  K 


The  binomial  and  negative  bi^nomial  expansions  are, 
respectively. 


(1-y)"  =  \    iTxSriyr  (-y)'  >  and 


Thus,  for  any  Integer 
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(1-xz) 


<^ 


1  =  0 


YiJpTTT 


(-1)  (xz)   ,  and 


-2( 


After  multiplying  these  quantities,  we  obtain 


(1-xz) 
(1-x) 


1-xz)'^   =   V   Y  r  i^i    (2(})+t-l)  icj)!     1  t+1 
l-x)2*   '  1^0  tio^^   l!(*-l)!t!(.2^-l)!  ^  ^ 

(3.3.7) 


Letting  S(z)  =  z"^"^"'"  ^^  ^^l^      and  denoting  by  C'',"the 

(l-x)2*  ^' 

coefficient  of  x  '  In,"  we  obtain 


c'^^^Siz)    ^    ^ 


I    (_n  n1       (j)  !  (2(j)+m-l)  ! m+1+1 

1  =  0^  "-^       1!  ((j)-l)!  (m+l-l)!(2c()-l)  !  ^ 


m+1 


(j)!  (2(})+m-l)  ! 


^io^  ^^   i!(4'-l)  !(m+l-l)!(2<|)-l)! 


((|)  <  m+1) 

m+1  +  1 

(cf)  >  m+1) 
(3.3.8) 


and, 


9^m+l„ 


r^^iuxi  ,  ,,,  9i  ^    Y  .  ,vl  (2({)+m-l)  !  (m+1+1) 

^9z^x   ^^-^^lz  =  l    (2cj)-l)!  J,^-^^   l!((|)-l)!(m+l-l)! 


where 


1  =  0 
'(})     If   (j)  <  m+1 

m+1   If   A  >  m+1 


(3.3.9) 


But, 
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S(z) 


z  =  l 


l9z       (i-x)2^ 


z=l 


-((f)+l) 

=  (1-x)       {  (m+l)-x(in+l+(j))} 


{(m+l)-x(m+l+(j))}   ^  ^j'l^?'  x^ 
i=0  ^- ^' 


so  that 


.m+1 


3z 


S(z) 


z  =  l 


.  ^  (m+l)  ((|)+m+l)  !  _  (m+l  +  ij))((|)+m)  ! 
(m+l)  Icj)!  m!(})  ! 


(3.3.10) 


=  0  . 
Equating  (3-3.9)  and  (3. 3. 10),  we  have 

r 

V  f    1  ^i  (2(|)+m-l)  !  (m+i-f-1)     n         x         •  /a   _i_->  i 
Z  (-1)   i!(c{)-i)!(m+l-i)!  -  =  0  '  6=min{(J,,m+l} 
1=0  vy   ;  ^       ;  (3.3.11) 

By  letting  m  =  k,  (j)  =  j-k  in  (3.  3. 11),  we  get  C  =  0 

in  (3.3.6)  demonstrating  that  Q      given  by  (3-3.1)  is 

indeed  the  inverse  of  Q.      The  matrix  Q~      is  now 

completely  specified;  any  element  may  be  obtained 

through  a  simple  arithmetic  calculation.   The  following 

recursive  method  of  generating  the  elements  of  Q 

may  also.be  useful: 

■  ^^"'^m,k  =  ^2-^)^^"')m-l,k-l--^-2)(^"')m-2,k-l  • 

(3.3.12) 


This  is  verified  by  obtaining  from  (3-3.1), 
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(«■').  1.1=  ./-^^"""^^"-^^'"'-     ,  and 

in-l,ic-i  2'""^(m-k)  !  (k-1)  !  (2k-m)  ! 

(Q-^)  =  (_i)"^-i<-l(2k-2)!-(^-i)!    ^ 

m-2,k-l  2™"^(m-k-l) !(k-l) ! (2k-m+l)! 


Now,  the  right-hand  side  of  (3.3.12)  is 

(-l)"^~^"(2k-2)  !  (m-1)  !     .m(2m-l)  ^  m(m-2)-. 
2">-2(m-k-l)!(k-l)!(2k-m)!   ^(m-k)    2k-m+l 

(-l)'^"^(2k)  !  (m+1)  !   ^(m-k)k(2k-m+l)     (m+l)(2k-l) 

o^/   ,^,1,/ol   ^1  ^  I  ""  2k(2k-l)(m+l)  ""  2  (m-k)  ( 2k-m+iy 
2  (m-k) ! k ! ( 2k-m+l ) ! 

_   (-l)"'"^(2k)!(m+l)!  ^  ^^-1) 


2'"(m-k)  !k!(2k-m+l)  !         ^'^ 

again  using  (3.3-1).   In  addition,  a  recursive  formula 
for  obtaining  the  elements  of  Q   is 

Using  (3.3.2)  we  have 


(2m-k)  ,    .  ^    (2m-k)  (m-1)  !  (2m-k-l)  !2^~-^ 

k(2m-l)  ^^V-l,k-l    k(2m-l)(2m-2)  !  (k-1)  !  (m-k)! 


m!  (2m-k)  !     2^(—  -) 


(2m) !k! (m-k) !    ^  m  2 


m!    (2m-k) !  pk 
(2m) !  k! (m-k) ! 


=  ^«)m,k 
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verifying  (3.3.13)  • 

3 . ^   Examples 

The  following  examples  will  Illustrate  the  devel- 
opments of  the  two  preceding  sections. 

1.   Let  it  be  required  to  find  the  distribution  of 
_  1  '    1  " 

1     ~  "o     o  o"  o   • 

The  characteristic  function  of  T  is 

-l^^el  2 

<l>n^iQ)    =    (e  {l+\hQ\)}  (using    (2.2.7) 

-|0|  ,  1         2 

=    e  (1  +  |e|   +^|9r) 

=  ^  e  (1  +  |e|)    +  I  e  (1  +  |6|  +^|e|^) 

=  }   ^3(6)  +  I  4.^(6)  . 
Thus,  Pr{T<t}  =  ^Pr{X  <t}  +  |pr{X  <t} 
and  the  density  of  T  is  given  by 

ii.   It  is  desired  to  find  the  distribution  of 

_  5      1     1 
T  -  j^X^  +  3-X3  +  j^X^  . 

The  characteristic  function  of  T  is 

(t>^(e)  =  (f)(^e;l)-(},(ie;3)-4>(^6;5) 

-   e    (i  +  i|9|)(i  +  i|e|  +^|e|2) 
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-|e| 

=  e     £'0  using  (3-2.4). 

Thus,  V    =    {    1      ^     ^     ^      0   0   O---) 
Using  Table  3-2,  we  get 


n-1^     -^iiririprTwiirir"" 


so  that 

Pr{T  <  t}  =  Y^VriX^   <  t}  +  p^Pr{X^  <  t}  +  ^Pr{X^  <  t} 
+  3^Pr{X^<t}  . 


111.   Consider  the  case  when  T  represents  the  average 

of  n  lid  pre-t  r.v.'s,  each  with  2m+l  d.f.   It  may  be 

st 
verified  that  for  n  =  2,  the  (1+1)    non-zero  element 

of  n  Is  given  by 


(Xn.     (2m+21)  !  (2m-21)  ! 
^1^    (m+1) ! (m-1) ! 


,  1=0,1,2, . . . ,m 


^"^"^    y  (m.    (2m+2j)!(2m-2j)! 

^to^  (m+j)!(m-j)!  (3.4.1) 


For  m=0, 1, 2, . . . , 7 ,  Table  3-3  gives  the  values  of  n. 

Iv.   It  will  be  convenient  for  future  reference  to 
calculate  the  jn-vectors  for  some  common  cases.   Letting 
n=2  In  (3^1.2),  we  have 


36 


Table  3.3 

THE  ELEMENTS  OF  n  FOR  THE  AVERAGE  OF  TWO  lid 
PRE-t  R.V.'s  WITH  2in  +  1  D.F. 


1      0 

0 

0 

0 

0 

0 

0  , 

.2500 

0 

0 

0 

0 

0 

0 

•  7500 

.0625 

0 

0 

0 

0 

0 

.2083 

.0156 

0 

0 

0 

0 

.7292 

.0656 

.0039 

0 

0 

0 

.1969 

.0201 

.0010 

0 

0 

.7219 

.0663 

.0060 

.0002 

0 

.1915 

.0222 

.0017 

.0001 

.7182 

.0665 
.1884 
.7160 

.0072 
.0234 
.0666 
.l86iJ 
•  71^5 

.0005 
.0023 
.0080 

.0241 
.0666 
.1850 
.7134 
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T  -  aX    +  (l-a)X    .  (3.^,1) 

1  2 

POT^  rv-9    ^32111   1      .1 

"  ~  T^'  5'  IT'  3'  2'  3'  IT'  5'  ^""'^  I^'  ^^^   ^°^ 

v^  =  1(2)9,  V2  =  v-j^(2)9,  n  was  computed  and  tabulated 
in  Table  3.  •4.   The  upper  percentage-points  of  the 
distribution  of  T  were  then  calculated.   For  (i>  =   .10, 
.05,  .025  and  .01,  Table  3-5  gives  values  of  A  (a,v^,V2) 
such  that 

Pr{T  >  X^{a,v^,v.^)}  =  <^    . 

Thus,  for  example, 

Pr{|x2  +|x     >  1.277}    =    .025    . 

3 . 5   The  Case  of  Even  Degrees  of  Freedom 

We  have  thus  far  restricted  attention  to  t  or 
pre-t  r.v.'s  having  odd  d.f.   This  has  been  necessitated 
by  the  form  of  ^^j^i^   in  (2.2.5),  which  gives  no  indica- 
tion of  the  form  of  (t>  ( 0  ;  v )  v/hen  v  is  even.   The  character- 
istic function  of  a  t  r.v.  vjith  even  d.f.  is  unknown. 

If  T  in  (3.1.2)  contains  a  X   for  v/hich  v.  is 
even,  we  may  put  bounds  on  its  distribution  using  our 


knowledge  of  the  distribution  of  T  when  v.  is  odd  for 
all  j.   To  do  this  we  need  some  properties  of  stochastic 
ordering. 

A  r.v.  X  is  said  to  be  stochastically  larger  than 
the  r.v.  Y  (X>Y)  iff  Pr{X  >  a}  >  Pr{Y  >  a}  for  all  a. 
VJe  may  define  an  "absolute  stochastic  ordering"  in  a 
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Table  3.4 


ELEMENTS  OF  n  WHEN  T =  aX   +  (l-a)X 

^1        ^ 


(V^.V^)   9/10      h/5       iJH  2/3       i/2       I/3       \/h  1/5       l/lO 


(1 

,3) 

.9000 

.8000 

.7500 

.6667 

.5000 

.3333 

.2500 

.2000 

.1000 

.1000 

.2000 

.2500 

.3333 

.5000 

.6667 

.7500 

.8000 

.9000 

(1 

.5) 

.9000 

.8000 

.7500 

.6667 

.5000 

.3333 

.2500 

.2000 

.1000 

.0900 

.1600 

.1875 

.2222 

.2500 

.2222 

.1875 

.1600 

.0900 

.0100 

.0400 

.0625 

.1111 

.2500 

.4444 

.5625 

.6400 

.8100 

(1 

,7) 

.9000 

.8000 

.7500 

.6667 

.5000 

.3333 

.2500 

.2000 

.1000 

.0882 

.1563 

.1781 

.2074 

.2250 

.1926 

.1594 

.1344 

.0738 

.0108 

.0384 

.0563 

.0889 

.1500 

.1778 

.1688 

.1536 

.0972 

.0010 

.0080 

.0156 

.0370 

.1250 

.2963 

.4219 

.5120 

.7290 

(1 

,9) 

.9000 

.8000 

.7500 

.6667 

.5000 

.3333 

.2500 

.2000 

.1000 

.0874 

.1509 

.1741 

.2010 

.2143 

.1799 

.1473 

.1234 

.0669 

.0112 

.0384 

.0552 

.0847 

.1339 

.1481 

.1356 

.1207 

.0729 

.0013 

.0091 

.0167 

.0353 

.0893 

.1411 

.1507 

.1463 

.1039 

.0001 

.0016 

.0039 

.0123 

.0625 

.1975' 

.3164 

.4096 

.6563 

(3 

3) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.7300 

.5200 

.4375 

..3333 

.2500 

.3333 

.4375 

.5200 

.7300 

.2700 

.4800 

.5625 

.6667 

.7500 

.6667 

.5625 

,4800 

.2700 

(3 

5) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.7290 

.3120 

.4219 

.2963 

.1250 

.0370 

.0156 

.0080 

.0010 

.2260 

.3280 

.3438 

.3333 

.2500 

.2222 

.2813 

.3520 

.5940 

.0450 

.1600 

.2344 

.3704 

.6250 

.7407 

.7031 

.6400 

.4050 

(3, 

7) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.7290 

.5120 

.4219 

.2963 

.1250 

.0370 

.0156 

.0080 

.0010 

.2186 

.3072 

.3164 

.2063 

.1875 

.0741 

.0352 

.0192 

.0027 

.0461 

.1360 

.1797 

.2364 

.2500 

.1975 

.2109 

.2560 

.4860 

.0063 

.0448 

.0820 

.1728 

.4385 

.6914 

.7383 

.7168 

.5103 

(3, 

9) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.7290 

.5120 

.4219 

.2963 

.1250 

.0370 

.0156 

.0080 

.0010 

.2157 

.2984 

.3051 

.2822 

.1741 

.0670 

.0314 

.0170 

.0026 

.0462 

.1317 

.1695 

.2116 

.2009 

.1058 

.0565 

.0329 

.0021 

.0091 

.0464 

.0771 

.1358 

.2188 

.1975 

.1846 

.2048 

.4084 

.0000 

.0115 

.0264 

.0741 

.2812 

.5926 

.7119 

.7373 

.5859 

(5, 

5) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.5905 

.3280 

.2383 

.1358 

.0625 

.1358 

.2383 

.3280 

.5905 

.3150 

.37  33 

.3516 

.2881 

.2083 

.2881 

.3516 

.3733 

.3150 

.0945 

.2987 

.4102 

.5761 

.7202 

.5761 

.4102 

.2987 

.0945 

Table  3.4  -  Continued 
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(v^.v^) 

9/1 0 

4/5 

3/^ 

2/3 

1/2 

1/3 

!/■* 

1/5 

1/10 

(5,7) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.5905 

.3277 

.2373 

.1317 

.0313 

.0041 

.0010 

.0003 

.0000 

.3012 

.3408 

.3115 

.2346 

.0938 

.0741 

.1377 

.2112 

.4784 

.0913 

.2240 

.2665 

.2881 

.2187 

.2305 

.3076 

.3584 

.3685 

.0170 

.1075 

.1846 

.3457 

.6562 

.6913 

.5537 

.4301 

.1531 

(5,9) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.5905 

.3277 

.2373 

.1317 

.0313 

.0041 

.0010 

.0003 

.0000 

.2952 

.3277 

,2066 

.2195 

.0781 

.0137 

.0036 

.0013 

.0000 

.0916 

.2102 

.2419 

.2455 

.1406 

.0604 

.0857 

.1382 

.3876 

.0200 

.1006 

.1516 

.2222 

.2344 

.1975 

.2571 

.3195 

.3958 

.0027 

.0338 

.0725 

.1811 

.5156 

.7243 

.6526 

.5407 

.2166 

(7,7) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.4783 

.2097 

.1335 

.0590 

.0156 

.0590 

.1335 

.2097 

.4783 

.3482 

.3172 

.2574 

.1575 

.0656 

.1575 

.2574 

.3172 

.3482 

.  1  398 

.2839 

.3045 

.2765 

.1969 

.2765 

.3045 

.2839 

.1309 

.0337 

.1892 

.3045 

.5070 

.7219 

.5070 

.3045 

.1892 

.0337 

(7,9) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0   . 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.4783 

.2097 

.1335 

.0585 

.0078 

.0005 

.0000 

.0000 

.0000 

.3396 

.3020 

.2419 

.1431 

.0313 

.0283 

.0755 

.1343 

.3874 

.1379 

.2585 

.2656 

.2189 

.0844 

.1021 

.1957 

.2710 

.3698 

.0380 

.1595 

.2175 

.2656 

.2062 

.2414 

.3046 

.3136 

.1865 

.0062 

.0703 

.1414 

.3139 

.6703 

.6277 

.4242 

.2811 

.0563 

(9,9) 

0 

0 

0 

0 

0 

0 

0 

0 

0  , 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.3874 

.1342 

.0751 

.0261 

.0038 

.0261 

.0751 

.1342 

.3874 

.3585 

.2554 

.1811 

.0859 

.0201 

.0859 

.1811 

.2554 

.3585 

.1816 

.2770 

.2540 

.1707 

.0663 

.1707 

.2540 

.2770 

.1816 

.0605 

.2129 

.2626 

.2690 

.1916 

.2690 

.2626  , 

.2129 

.0605 

.0120 

.1205 

.2272 

.4484 

.7182 

.4484 

.2272 

.1205 

.0120 

Table  3.5 

VALUES  OF  A  SUCH  THAT  PR{aX   +(l-a)X   >A}  =  d) 

^1    .    ^2 
((v.,v„)  given  in  parentheses) 
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.10 


,05 


.025 


,01 


a 

(1 

.3) 

.900 

2.774 

5.684 

11.437 

28.640 

.800 

2.478 

5.060 

10.170 

25.459 

.750 

2.335 

4.751 

9.538 

22.869 

.667 

2.106 

4.242 

8.488 

21.222 

.500 

1.695 

3.267 

6.416 

15.936 

.333 

1.366 

2.413 

4.444 

10.701 

.250 

1.234 

2.062 

3.561 

8.154 

.200 

1.165 

1.881 

3.093 

6.689 

.100 

1.044 

1.583 
(1 

2.342 
5) 

4.157 

.900 

2.771 

5.683 

11.436 

28.640 

.800 

2.467 

5.054 

10.166 

25.458 

.750 

2.317 

4.740 

9.532 

23.867 

.667 

2.069 

4.218 

8.476 

21.217 

.500 

1.593 

3.185 

6.367 

15.916 

.333 

1.178 

2.195 

4.291 

10.728 

.250 

1.006 

1.743 

3.251 

7.980 

.200 

0.918 

1.508 

2.665 

6.402 

.100 

0.772 

1.142 

1.689 

3.339 

10 


,05 


,025 


01 


(1,7) 
,900  2.771        5.683      11.436      28.640 

,800  2.465        5.053      10.166      25.458 


750 

2 

313 

4 

738 

9.531 

23 

867 

667 

2 

062 

4 

215 

8.474 

21 

217 

500 

1 

570 

3 

173 

6.361 

15 

914 

333 

1 

114 

2 

149 

4.256 

10 

615 

250 

0 

919 

1 

659 

3.201 

7 

891 

200 

0 

819 

1 

394 

2.597 

6 

385 

100 


,900 
,800 


0.655        0.971 


1.485 


3.241 


(1,9) 
2.771        5.683      11.436      28.640 
2.464        5.052      10.166      25.458 


750 

2.314 

4 

747 

9.569 

24 

108 

667 

2.059 

4 

213 

8.473 

21 

217 

500 

1.561 

3 

168 

6.359 

15 

913 

333 

1.088 

2 

139 

4.267 

10 

719 

250 

0.877 

1 

633 

3.202 

7 

966 

200 

0.766 

1. 

347 

2.579 

6 

379 

100 


0.588   0.880   1.399   3.220 


(3 

,3) 

0.861 

1.234 

1.664 

2.369 

0.798 

1.135 

1.519 

2.145 

0.775 

1.095 

1.460 

2.050 

0.746 

1.047 

1.385 

1.924 

0.723 

1.009 

1.324 

1.821 

0.746 

1.047 

1.385 

1.924 

0.775 

1.095 

1.460 

2.050 

0.798 

1.135 

1.519 

2.145 

0.861 

1.234 

1.664 

2.369 

(3 

.5) 

0.855 

1.227 

1.657 

2.362 

0.773 

1.103 

1.484 

2.109 

0.736 

1.046 

1.401 

1.979 

0.683 

0.961 

1.277 

1.790 

0.606 

0.835 

1.081 

1.463 

0.573 

0.778 

0.989 

1.295 

0.573 

0.777 

0.983 

1.278 

0.579 

0.785 

0.996 

1.293 

0.607 

0.827 

1.052 

1.373 

(3, 

7) 

0.853 

1.225 

1.-656 

2.361 

0.766 

1.097 

1.478 

2.104 

0.726 

1.035 

1.392 

1.978 

0.663 

0.939 

1.254 

1.771 

0.562 

0.777 

1.012 

1.388 

0.501 

0.676 

0.853 

1.111 

0.487 

0.653 

0.816 

1.042 

0.485 

0.648 

0.808 

1.024 

0.497 

0.664 

0.827 

1.046 

(3, 

9) 

0.853 

1.224 

1.655 

2.361 

0.764 

1.094 

1.476 

2.102 

0.721 

1.031 

1.388 

1.974 

0.654 

0.929 

1.245 

1.764 

0.540 

0.749 

0.983 

1.362 

0.461 

0.622 

0.788 

1.035 

0.439 

0.585 

0.729 

0.928 

0.431 

0.572 

0.708 

0.891 

0.432 

0.572 

0.705 

0.878 

Table  3.5  -  Continued 


^1 


10 


,05 


025 


01 


a 

(5 

.5) 

900 

0.599 

0.817 

1.040 

1 

.360 

800 

0.550 

0.746 

0.946 

1 

.229 

750 

0.531 

0.717 

0.905 

1 

.169 

667 

0.506 

0.680 

0.854 

1 

.095 

500 

0.487 

0.651 

0.811 

1 

029 

333 

0.506 

0.680 

0.854 

1 

095 

250 

0.531 

0.717 

0.905 

1 

169 

200 

0.550 

0.746 

0.946 

1 

229 

100 

0.599 

0.817 
(5 

1.040 
7) 

1 

360 

900 

0.597 

0.815 

1.038 

1 

358 

800 

0.542 

0.736 

0.935 

1 

218 

750 

0.518 

0.701 

0.888 

1 

153 

667 

0.483 

0.650 

0.817 

1 

051 

500 

0.439 

0.584 

0.724 

0. 

912 

333 

0.432 

0.572 

0.706 

0. 

882 

250 

0.443 

0.587 

0.727 

0. 

912 

200 

0.454 

0.604 

0.749 

0. 

943 

10 


05 


025 


.01 


100 


0.487   0.652   0.812 


1.028 


(7 

,7) 

0.485 

0.649 

0.809 

1.024 

0.444 

0.592 

0.735 

0.927 

0.428 

0.569 

0.705 

0.887 

0.406 

0.538 

0.663 

0.827 

0.389 

0.513 

0.629 

0.777 

0,406 

0.538 

0.663 

0.827 

0.428 

0.569 

0.705 

0.887 

0.444 

0.592 

0.735 

0.927 

0.485 

0.649 

0.809 

1.024 

(7 

9) 

0.484 

0.648 

0.807 

1.023 

0.440 

0.587 

0.729 

0.921 

0.421 

0.560 

0.695 

0.876 

0.393 

0.521 

0.644 

0.805 

0.362 

0.475 

0.581 

0.716 

0.362 

0.476 

0.582 

0.715 

0.375 

0.494 

0.606 

0.749 

0.387 

0.511 

0.628 

0.778 

0.419 

0.555 

0.684 

0.852 

.900 
.800 
.750 
.667 
.500 
.333 
.250 
.200 
.100 


(5,9) 
0.596        0.814        1.037        1.375 
0.538        0.732        0.931        1.214 


0.511 
0.471 
0.414 
0.389 
0.392 
0.398 


0.694 
0.636 
0.550 
0.513 
0.516 
0.524 


0.881 
0.802 
0.682 
0.630 
0.633 
0.644 


1.149 
1.038 
0.862 
0.779 
0.782 
0.798 


0.422        0.558        0.688        0.856 


42 

similar  manner: 

Definition :   The  r.v.  X  is  said  to  be  absolutely 
stochastically  larger  than  the  r.v,  Y  (denoted  X^Y) 
iff  Pr{|x|  > a}  >  Pr{|Y|  >  a}  for  all  a  (andstrict 
inequality  for  at  least  one  a). 

Notice  that  x'^  Y  iff  |x|  f  |Y|  . 
Lemma 

Let  X  and  Y  be  independent  random  variables  with 
x'^Y.   If  Z  is  independent  of  X  and  Y  then  X  +  z'^Y  +  Z. 
proof: 
x'^  Y      =>        \X\  %   |Y| 

=»  Pr{|X|  >  u}  >  Pr{|Y|  >u},  for  all  u 
=>  Pr{|x|  <  u}  <  PrlJY]  <u},  for  all  u 
=i>       Pr{-u<X<u}  <  Pr{-u<Y<u},       (3-5.1) 


for  all  u, 


Let    Z    have    cdf   F    .       Then, 


Pr{-w  <  X  +  Z  <  w}    = 


Pr{-w-C  <  X  <  w-C}dF    (^ 


and. 


Pr{-w  <  Y  +  Z  <  w}    = 


Pr{-w-C  <  Y  <  w  C}dF    (K) 


From    (3.5.1)    we   have 

Pr{-w-C  <  X  <  w-O    ^  Pr{-w-5  <  Y  <  w-U,    for   all   w,^ 

=>     Pr{-w<X  +  Z<w}    <  Pr{-w  <  Y  +  Z  <  w}  ,    for   all   w 

=>lX  +  Z|>^|Y  +  Zl  =>       X  +  z'^Y  +  Z  □ 


^3 

THEOREM    3-1 

Let    X.'f'Y.,    1  =  1,2,    where    X,,    X^ ,    Y,    and  Y„    are 
11  '     '  1'       2 '       1  2 

independent   r.v.'s.  Then   X.  +Xp>Y,  +Yp. 

proof: 

X^  +  X^'^'y^  +  X^  by    Lemma 

Y^  +  X^  '^'  Y^  +  Yp  by    Lemma 

X^ +X2'^' Y^  +  Yp  by    transitivity  U 

Is! 
From    §2.1    we    see    that    X       >  X        whenever   v.,  <  v„  . 

v^    Vp  12 

Therefore,  from  the  theorem,  the  distribution  of 

n 
T=  y  a.X   wherein  some  of  the  v.'s  are  even,  may  be 

•   -,1V.  1  }  J 

j=l  ^   J  ^   n 

bounded  by  the  distribution  of  T,  =  T  a.X  i  and 

T^  =   )  a . X  "  where 

V.  =  V.  =  V.  if  V.  is  odd,  and 

J    J    J      J 

v.=v.-l,   v.=v.+l  if  V.  is  even. 
For  example  should  one  desire  the  distribution  of 

T  =  aX^  +  (1  -  a)X^  , 

he  may  find  the  distributions  of  T-,  =  aX^  +  (1  -  a)X|- 

i     j         D 

and  Tp  =  aX  +  (1  -  a)Xt-  and,  noting  that  X  '^'xJf^X  , 

Isl  Isl 
Theorem  3.1  implies  that  T  '>T>Tp,  or  alternatively, 

Pr{  |t^I  >  x}  >  Pr{  |aX^  +  (1  -  a)Xt-|  >  x}  >  Pr{  JTpl  >  x}, 

for  all  x. 

For  one-sided  approximations,  vje  may  remove  the 

absolute  value  signs  since  T,  ,  Tp  and  T  are  each 

symmetrically  distributed  about  0. 


CHAPTER  4 
THE  POWER  OF  THE  T-STATISTIC 


4  . 1   Introducti( 


Let  u  (nxl)  be  a  vector  of  independently  distri- 
buted random  variables  with  knovm  distribution  and 
suppose  X  and  u  are  related  by 

X  =  u  +  dl_  , 

where  1  represents  the  (nxl)  vector  of  unit  elements. 

Given  X,  suppose  it  is  .required  to  distinguish  between 

HqI  d=0    and 

K,  :  d>0 
± 

subject  to  Prdncorrectly  rejecting  H  }  =  (}> .   VJe  consider 

n 
tests  of  the  type  a'x=  I   a.x..   Given  d,  let  T,(a) 

1=1  ^  ^  ^ 

denote  a  r.v.  with  the  same  distribution  as  a'x.   That 

is, 

n  n       n 

T  (a)  ~  I   ax   =  d  ^  a   +  I   a  u  . 
1=1  ^  ^    1=1  ^    1=1  ^  ^ 

Defining  A,  (a)  to  be  the  critical  value  for  rejection 
in  a  size  <p    test,  we  have 

Pr{TQ(a)  >  X^(a)}  =   ^  . 
The  power  function  of  the  test  is 

Pr{T  (a)  >  A, (a)}  .  (4.1.1) 


■d 
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^5 


The  problem  now  confronting  us  is  that  of  choosing 
a  to  maximize  Cl.l.l),   If  the  u.'s  are  symmetrically 
distributed  about  0,  the  distribution  of  T„(a)  is 
invariant  over  sign  changes  of  the^  elements  of  a.  Vie 
may  assume  therefore  that  a.  >  0,  i=l,2,...,n.   In 
addition,  for  any  positive  constant  k,  we   have 


Pr{TQ(Ka)  >  A^(K:a)}  =  Pr  {T^  (a)  >  A  (a ) }  = 


But 


and 


so  that 


and 


TQ(Ka)  =  KTpCa) 


T^(Ka)  =  KT^(a) 


A,(Ka)  =  KA,(a) 


Pr{T^(Ka)  >  A  (Ka)}  =  Pr{T^(a )  > ^A^ ( Ka) } 

-  Pr{T^(a)  >  A^(a)}. 

n  n 

We  may  therefore  let  I   a.  = 1,  for  if  not,  let  k  =  1/  I  a.  . 

i=l  ^  1  ^ 

Our  attention  is  restricted  to  the  case 

n 
a.  >0,  i  =  l,2,...,n  and  "[   a.  =1. 

i=l 


If  n = 2  and  u.  =X   ,  the  statistic  T,(a)  becomes 
1    V.  '  d  — 

the  test  criterion  of  §1.2.   As  a  special  case  of  the 
above,  therefore,  is  the  choice  of  weights  assigned  to 
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the  t-variable  needed  to  maximize  the  power  of  the 
test  used  to  solve  the  Behrens-Pisher  problem.   This 
chapter  gives  some  general  theorem.s  on  comparing  the 
power  of  tv;o  tests.   Some  examples  are  given,  along 
with  some  counterexamples  to  intuitive  notions.   It 
is  shown  that  the  power  of  the  pre-t  r.v.  increases 
with  d.f.  and  that  there  exists  no  distinct  set 
{a.},  i=l,2,...,n,  which  yields  maximum  power  of  T,(a) 
for  all  size  tests  when  the  d.f.  differ. 

4 .2   General  Results 

The. theorems  which  follow  give  a  general  method 
of  comparing  the  power  functions  of  tv;o  tests  about  a 
hypothesis  of  a  location  parameter  and  will  be  useful 
in  the  ensuing  Investigation. 

THEOREM  4.1 

Let  X  and  Y  be  absolutely  continuous  random  variables 
v/ith  distribution  functions  P  and  G  respectively,  and  let 
a  and  6  satisfy  P(a)  =G(3).   If  the  inverse  function  of  F, 
say  P  ,  exists,  then  the  following  conditions  are 
equivalent : 

(i)   P(a-d)  >  G(3-d)      V  d  >  0 
(li)   P  G(x+d)  >FG(x)+d    Vx,  d>0 

(iii)     dx    "  -^  ^   ^   where  P  G  is  dif f erentiable 

I  •     \       Sot  ^  -, 

(IV)   ^>1  .       ^ 
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proof: 
(1)    iff    (ii): 


F(a  -  d)    >   G(3  -  d),       Vd  >  0 


a-d>FG(3-d),       Vd>0 


a   =   F   G(6)    >    F   G(3-d)+d,  Vd>0 


F   G(x  +  d)    >    F   G(x)+d,  Vx,d>0    (6  =  x  +  d) 


(11)  Iff   (iii; 


F   G(x  +  d)    -   f'g(x)    >    d, 


F*G(x  +  d)    -    F*G(x) 


V  X  ,  d  >  0 


>    1,  V  X  ,  d  >  0 


dF    G(x) 
dx 


>    1,  V  X 


Conversely, 


(111) 


x+d       ^  rX+d 

dF   G(y)    >  dy 

X  •'x 


F'G(x+d)-F   G(x)    >    d  Vx,d>0, 


(111)    Iff    (Iv): 


a   =    F   G(3) 


da   _    dF    G(3) 
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d3 


>    1 


THEOREM  4.2 

Let  X  and  Y  be  absolutely  continuous  r.v.'s  with 
cdf's  F  and  G  respectively..  Suppose  a  and  g  satisfy 
F(a)  =  G(3).   Then,  Y  +  d  is  a  more  powerful  test  of  size 
4)  than  X  +  d  for  testing 

H^:  d=0 
vs.   H.  :  d  >  0 
iff  rr^  >  1 .      V/e  shall  indicate  that  this  property  holds 

0  P 

by  writing  F<G  (read  "F  is  less  in  power  than  G"). 
(P) 

proof: 

|£>1   <j=^   F(a-d)  >  G(3-d),     V  d>0 

<=?>   •Pr{X<a-d}  >  Pr{Y<B-d},    V  d>0 

<=^   Pr{X  +  d>a}  >  Pr{Y  +  d>3},    V  d>0 

<?=^   Y  +  d  is  more  powerful  than  X  +  d  in  a 
test  of  Hp  vs .  H-.  . 

If,  in  addition  to  the  hypotheses  of  Theorem  '^.l, 

there  exists  a  c  such  that  F(c)=G(c),  |c|  <  «>,  and  if 

any  of  (i)  thru  (iv)  in  Theorem  4.1  holds,  we  have 

a>6  (if  &  >  c) .      This  is  seen  by  writing 

FG(x  +  d)>FG(x)+d   =^  FG(c+d)>c+d,   Vd>0 


so  that 


a  =  F*G(3)  >  3 


The  next  result  gives  a  relationship  between  the 

power  of  two  tests  and  their  moments,  giving  a  necessary 

condition  for  F<G. 
(P) 
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THEOREM  i|  .  3 


Let  X  and  Y  be  absolutely  continuous  r.v.'s  distrib- 
uted symmetrically  about  0,  with  cdf's  F  and  G  respec- 
tively.  Suppose  that  F  Is  strictly  Increasing  so  the 
Inverse  function,  F  ,  exists  and  Is  uniquely  defined. 

If  F  <  G  then  E(Y^^)  <  E(X^^^),  k  =  1,2,3,  ••  •  providing  these 

(p) 
expectations  exist. 


proof: 


F  <  G 
(P) 


da 
96 


>  1 


Theorem  i| .  2 


X  =  0 


«  s 

F  G(x  +  d)>F  G(x)+d,  Vd>0,  x   Theorem  ^.1 

F  G(d)  >  F  G(0)  +d  =  d 


or     G(x)>F(x)      Vx>0 

=i>  1  -  G(x)  <  1  -  F(x)     Vx  >  0 


ilk 


•  oo  .  oo 

x2^-l[l  _G(x)]dx  <  4k   x^^~^[l  -F(x)]d: 


x^^^dG  <  2 


X   dF 


["'x^^dG  . 


0 
2k 


dF 


(integrate  by  parts) 


or,  ^(Y^^)  <  E{J^^) 


D 


4 • 3   Applications:   Some  Examples  and  Counterexamples 

Let  V/  be  a  r.v.  with  cdf  H(')  and  Z  a  constant  3 

0<Z<1.   If  H^C-)  denotes  the  cdf  of  V//Z,  then  H„   <  H„ 

^  ^l(p)  ^2 

whenever  Z^  <  Zp .   This  is  clear  upon  letting  H^  (a)  = 

9a   "2  ^ 


H„  (6),  so  that  aZ,  =  3Z„  or 


2 


>  1 


Now  invoke 


^1 
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Theorem  4.2.   This  situation  occurs  whenever  two 
random  variables  have  distribution  differing  only  in 

scale  parameters.   In  particular,  H„  <  H. 

^  (p) 
Suppose  now  that  Z  is  a  r.v.,  0<Z<1  with  prob- 
ability 1,  and  suppose  it  is  still  true  that 

H^  ^  H.  .         (4.3.1) 

(P) 

Of  primary  importance  in  the  investigation  into  the 
power  of  a  linear  combination  of  pre-t  r.v.'s  is  whether 
the  power  of  a  pre-t  r.v.  increases  with  d.f.   Let 


W=     Y 


(Xi  +  Xo) 


have  the  cdf  H  and 


Z  = 


<   2  ,    2y 


where  Y  is  a  standard  normal  r.v.  and  x-j  1=1,2  are 
independent  chi-square  r.v.'s  on  v.  d.f. 

W  then  has  the  pre-t  distribution  on  v  +v   d.f. 

and  is  independent  of  Z. 

Y 
W/Z  =  -—  has  the  pre-t  distribution  on  v^  d.f. 
Xi  1 

We  therefore  would  reach  the  conclusion  that  the  power 
of  the  pre-t  r.v.  will  increase  with  d.f.  if  (4.3.1) 
held. 
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It  Is  not  true  in  general,  however,  that  (4.3.1) 
holds.   The  following  is  a  counterexample. 

Let 

f    .4  ,  0  <  |xl   <  1 
h(x)    =    j    .1  ,  1 <  jx|  <  2      and      Z   = 
I       0  ,  elsewhere 


h     wp  h 

1      wp  \ 


Let   a,    3   be   3    H„(a)  =  H(3)  =  .95 


Then,    a  =  2. 0, 


1 . 5    since 


H^Ca)    =    E^viM  <  2Z} 


=    JgH(l)  ^\^{2)_ 

=    .95  =    H(1.5)    =    H(B) 


Taking   d  =  . 2    (say) , 


H^(a-d)    =    igH(l  -|)  +J2H(2  -d) 

=    i2H(.9)  +igH(1.8)    =    .92 


H(3  -d)    =   H(1.3)    =    .93 


so   that   H2(a  -d)  <H(3  -d),    or  ||  <  1   which  means    that 

H^  ^     H.     . 
(P) 

It  is  of  interest  to  investigate  the  conditions  for 

which  it  is  true  that  H„  <   H.   Letting  a,  6  be  3 

H2(a)  =  H(3),  we  have 


£'„H(aZ)  =  H(3) 


ff2Zh(aZ)da  =  h(3)d3    or 


1^  =    h(3) 
93    £'7Zh(aZ) 
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By  Theorem  4.2,  W       <     H  iff  h(3)  >S„Zh(aZ),   It  should 

(p)  ^ 

be  noted  that  H( B )  -  fi'^HCaZ )  < H(a)  so  that  B  < a.   The 

following  are  examples  for  which  H„  <  H. 

^  (P) 


(i)      Let   Z~Bln(l,p)    and   h(')    nonincreasi 


ng. 


H^Ca)    =    H(3) 


9a   _        h(6) 


h(a)  <  h(3) 
h(6)     ^    h(3) 


93        ff^ZhXaZT        ph(a)    ""   h(a) 


>    1 


(ii)      Let   h(5)  =  Ae    ^^,      ^  ,\  >  0 


H(3)    =    E^}]{aZ) 


^         -A3    _    T       p   ^-^aZ  -A3         „      -AaZ 

1  -  e  -    1  -  E„e  or,    e  =    E„e 


h(3)  =  Ae"^^  =  Aff^e"^^^  >  ff^AZe"^"^  =  ff^^hCaZ) 


3a   _         h(3) 


33         ff^^h^"^) 


>    1 


(iii)      Let    h(*)    be   non-increasing   and   v(0  =  ?h(?) 

concave 
.■      .  i.e.  ,     V6  e  (0,1)    and   any   ?^  ,  ^2 

6  v(q)  +  (1  -  6)  v(52^  ^    v(9C^  +  (1  -  6)^2)   . 

By   Jensen's    inequality,    ffv(aZ )  < v(aff(Z ) ) 
=   vCau)    where    \i  -  E(Z) . 

ff^ci^hCaZ)    <    ayh(ay)  '  _ 


ff^ZhCaZ)    <    yh(ay) 


E^maZ)    <    H(ay 
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But  H(3)  -  E^E(aZ)  <  H(ap) 


ay 


so  that 


Now. 


h(3)  >  h(ay) 


^   =    h(e)    >   h(g)   ^  h(g) 
9  3   ff^zFT^^  -  yh(ay)  "  h(ay) 


>  1 


An  example  of  such  a  distribution  is 
h(U  ^  ^^^  ~p^^  ,  0  <  5  <  K.   The  function  ?h(C)  is 
concave  since  — ^  ^  =  ^  <  0. 

(iv)   Let  h(0  =  pC^~^  ,  0  <  C  ^  1.   The  cdf  is 


H(0  = 


0  ,   5  <  0 
?^  ,  0  <  C  <  1 

1  ,   C  >  1 


'  ^z 

Zh 

aZ)  = 

|-min(l,  1/a 

0 

|-min(l,l 

a 
0 

) 
/ 

where 

P 

is  the 

cdf  of  Z. 

But 

H(3)  = 

ff^HCaZ) 

Zp(aZ)P  ^dP 


(aZ)^'dP 


'mind, 1/a)     . 
6^  -  I  (aZ)PdP  +Pr{Z  >  1/a} 


5^ 


min(l,l/a) 


(aZ)PdP 


so 


ff^Zh(aZ)    =   ^ 
z  a 


min(l, 1/a) 


(aZ)^dP 
..^bP   -<   §bP   =    pbP-I   =   h(B) 


or. 


3a 


>    1. 
op 

(v)      Let    h(C)  =  p(l  -  0^~^  ,  0  <  5  <  1  ,  p  >  1.      The   cdf   i; 

0  ,  ?<  0 
H(0    =    ^    1  -  (1  -  ?)P  ,        0  <  ?  <  1 

1  ,  C^  1 


H(3)    =    ^2^^"^) 


1  -  (1  -  3)' 


(1-  3)^  = 


mln(l,l/a) 


[1  -  (1  -  aZ)^]dP  +  Pr{Z  >  1/a} 


0 


mln(l,l/a) 


(1  -  aZ)^dP 


Now, 


E^Zh{aZ', 


■mind,  1/a) 


pZ(l  -  aZ)P~-^dP 


but 


min(l, 1/a) 


Z(l  -  aZ)P"^dP 
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min(l,l/a) 


Z^dP 


mln(l,  1/a)      p-1  p-1 

[(1  -  aZ)   ]   .dP 


P    p-1 


0 


by  Holder's  Inequality, 
Thus , 


£■  Zh(aZ) 


s  P 


^  P 


=  p 


min(l, 1/a) 


1 


zPdP 


min(l, 1/a) 


[(1  -  aZ)]^dP 


P-1 
P 


mln(l, 1/a) 


(1  -  aZ)PdP 


P-1 
P 


P-1 


(1  -0)^1  ^   =  p(l-6)^-"  =  h(6) 


P-1 


aa 


>  1 


4 .H      The  Power  of  the  Pre-t  R.V. 

Let  t^  have  the  Student-t  distribution  on  v  d.f. 
and  define  X^  =  t^//v.   For  v  -  2j  + 1,  j  =  0,1,2, ... , 
denote  the  cdf  of  X^  by  H  (^)  -  PrCX^  .^-^  <  U  ,  and  the 
pdf  by  h.  (O  =  tfH.(5)  .   Let  B   be  9  H  .  (  B. )  =  a.   We  wish 
to  demonstrate  that  H,   <  H  whenever  k  <  m.   It  will  be 

(p) 

sufficient  to  show  that 
93, 


k 


9  3. 


>  1 


Va 


and  then  use  Theorem  ^4 . 2 


56 


Let    k<m.      We   have,    from   Property    3   of    §2.4 


uja  =  H,  (u+c  I 


m-l   h. (C) 


j=k      '^ 


i^.H.l) 


so   that. 


m-l  dCh.C?) 


(i^.4.2) 


With 


•^1^^^    "         21  +1    1   '   ^1  +  5^)     ^"^^^^  >     Property    1    of    §2 


2     '2 


so   that. 


dCh    (O  2r(J  +l)h,(U 

il =    h,.(a ^ .  (4.4.3) 

a^  J  (1  +  ?    ) 


Notice   that 


dCh    (C)  p  P 

^1 >0        iff        (1  +  5)    s   2r(j  +1) 


iff        ^      < 


(2J  +1) 


(4.4.4) 


Now, 


H   (3    )   =  H,  (6,  ) 
mm  k      k 


m-l         ^        d?hj(C) 


k        mm' 


'^m~\^ 


W 


-,(^.^.5) 
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which  results  by  substituting  (H.H.2).      It  follows 
from  (2.1.2)  that  3j^<3^,  so  that  \(  3^)  ^  hj^(  g  )  . 
From  (^.^.5), 


>  1 


If 


>  0  ,   j  =k,k+l,. 


,m-l, 


By  (^.it.lj)  this  will  occur  Iff  |3„  < 
If  .5  s  a  <  .75,  then 


1 


/2j  +1 


H^(B^)  =  a  <  .75  =  Hq(1)  =  Pr{t^  <  1} 


<    Pr{t 


2m+l 


1}  =  H 


Therefore, 


1 


/2m  +  1 


/2m  +  1    /2j  +  1 
and  we  have 

!!k . 

9B™  ^ 


j  =  k,k+l, . . . ,m-l  , 


This  appears  to  be  true  also  when  o>   .75.   For  a  = 
.75(.05).95(.01).99  and  k  -  0, 1, 2 , . . . , 21  the  values  of 
B^  have  been  computed  and  tabulated  in  Table  4,1. 
Beneath  each  value  of  ^^   In  the  table  appears  the  value 
of  \(B^).   It  Is  seen  that  as  k  increases,  h,  (3,  ) 
Increases  also  for  each  a.   Thus, 
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Table   4.1 


VALUES   OF    B^(o)    FOLLOWED    BY   h    (3,  (a)) 


a 
'75  .80  .85  .90  .95  .96  .97  .98  .99 


0  1.000  1.376  1.962  2.800  2.800  2.800  2.800  2.800  2.800 
0.159  0.110  0.066  0.036  0.036  0.036  0.036  0.036  0.036 

1  0.441  0.565  0.722  0.946  1.359  1.504  1.704  2.010  2.621 
0.446  0.366  0.275  0.177  0.079  0.060  0.042  0.025  0.010 

2  0.325  0.411  0.517  0.660  0.901  0.980  1.083  1.233  0.505 
0.628  0.531  0.417  0.287  0.143  0.113  0.083  0.053  0.024 

3  0.269  0.339  0.423  0.535  0.716  0.773  0.847  0.951  1.133 
0.770  0.660  0.527  0.373  0.195  0.156  0.117  0.077  0.037 

4  0.234  0.294  0.367  0.461  0.611  0.657  0.717  0.799  0.940 
0.891  0.768  0.620  0.444  0.238  0.193  0.146  0.098  0.049 

5  0.210  0.264  0.328  0.411  0.542  0.581  0.632  0.702  0.819 
0.998  0.864  0.700  0.507  0.276  0.225  0.172  0.117  0.059 

6  0.193  0.241  0.300  0.374  0.491  0.527  0.571  0.633  0.735 
1.093  0.950  0.773  0.563  0.311  0.254  0.195  0.134  0.069 

7  0.179  0.224  0.277  0.346  0.453  0.485  0.525  0.581  0.672 
1.182  1.029  0.841  0.614  0.342  0.281  0.217  0.149  0.077 

8  0.167  0.210  0.259  0.323  0.422  0.452  0.489  0.540  0.622 
1.264  1.101  0.902  0.663  0.371  0.305  0.236  0.163  0.085 

9  0.158  0.197  0.244  0.305  0.397  0.424  0.459  0.506  0.583 
1.343  1.171  0.961  0.706  0.397  0.328  0.254  0.176  0.092 

10          0.150  0.187  0.232  0.289  0.376  0.402  0.434  0.478  0.549 

1.415  1.236  1,015  0.749  0.423  0.349  0.271  0.188  0.099 


Table  4.1  -  Continued 
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75 


,80 


.85 


90 


a 
.95 


96 


,97 


,99 


11 


0.143     0.179      0.221      0.275      0.357      0.382     0.413      0.454      0.521 
1.486      1.296      1.068      0.789      0.447      0.369      0.287      0.200     0.106 


12 


0.137      0.171      0.212      0.263      0.342      0.365      0.394     0.433     0.497 
1.552      1.356      1.118      0.826      0.470     0.388      0.303     0.211      0.112 


13  0.132  0.164  0.203  0.253  0.328  0.350  0.378  0.415  0.476 
1.615  1.414  1.164  0.863  0.493  0.407  0.317  0.221  0.118 

14  0.127  0.158  0.196  0.244  0.316  0.337  0.363  0.399  0.457 
1.678  1.469  1.209  0.897  0.512  0.425  0.331  0.232  0.124 


15 


16 


17 


18 


19 


20 


21 


0.123  0.153  0.189 

1.735  1.520  1.253 

0.119  0.148  0.183 

1.794  1.570  1.297 


0.235      0.305      0.325      0.351      0.385      0.441 
0.931      0.532     0.442     0.345     0.241      0.129 


0.228 
0.965 


0.295     0.314     0.339      0.372     0.426 
0.552      0.458      0.358      0.250      0.134 


0.115  0.144  0.178  0.221  0.285  0.305  0.328  0.361  0.412 

1.847  1.620  1.339  0.995  0.572  0.474  0.371  0.259  0.139 

0.112  0.140  0.173  0.214  0.277  0.296  0.319  0.350  0.400 

1.901  1.667  1.378  1.027  0.589  0.490  0.383  0.269  0.144 

0.109  0.136  0.168  0.209  0.270  0.288  0.310  0.340  0.388 

1.955  1.713  1.418  1.056  0.608  0.504  0.395  0.277  0.149 

0.106  0.133  0.164  0.203  0.263  0.280  0.302  0.331  0.378 

2.006  1.761  1.453  1.084  0.623  0.519  0.406  0.285  0.154 

0.104  0.130  0.160  0.199  0.257  0.273  0.295  0.323  0.368 

2.055  1.804  1.492  1.110  0.640  0.532  0.416  0.293  0.158 


60. 


>   1 


for  m  >  k, 


(4.4.6) 


Since  this  holds  for  all  o,    we  may  v/rite 


H,  <  H 


for  m  >  k, 


(4.4.7) 


It  is  Interesting  to  note  that  a  stronger  condition 
is  actually  true,  namely  the  function  H,  H  (x)  is  a 
convex  function  of  x.   That  this  implies  (4.4.6)  follows 
from  Theorem  4.1  since,  if  H,  "h   is  convex  then 


dH,  H  (x) 
k  m 

dx 


V  x. 


(providing 


dH,  H  (x) 
k  m 

dx 


,^n  >  1).   For  k  =  0,1,2,. ..,5  , 


^k  ^k+1*-^^  ^^^  been  computed  for  several  values  of  x, 


The  values 


^k  "k  +  l^^^^)-Hk  \+l^^) 
h 


are  seen  to  increase  with  x.   These  values  have  been 
tabulated  in  Table  4.2.   In  addition,  for  small  values 
of  X  the  tables  indicate  H,  H,  ^-,  (x)  >x  or, 

K    K  +  1 


\     "k  +  l^^^-«k  Hk-fl^°) 


h 


^  1  , 


supporting  the  belief  that 


dx 


x  =  0 


>  1 


Table  4.2 

VALUES  OF  H,  H,  ^,  (x) 
k   k+1  ' 

k  =  0  ^ 

\  "k+1^^^  FIRST  DIFFERENCE 

0.0004 
0.0997 
0.1013 
0.1031 
0.1061 
0.1104 
0.1152 
0.1216 
0.1287 
0.1371 

0.1466 

0.1570 

0.1690 

0.1818 

0.1957 

0.2109 

0.2274  ■ 

0.2447 

0.2639 

0.2834 

0.3045 

0.3265 

0.3503  ■ 

0.3743 

0.4004 

0.4271 

0.4551 

0.4844 

0.5145 
0.5457 

0.5788 
0.6121 
0.6470 
0.6829 
0.7198 
0.7585 
0.7977 
0.8381 
0.8800 
1.95  14.7235  0.9221 


0.00 

0.0004 

0.05 

0.1001 

0.10 

0.2014 

0.15 

0.3046 

0.20 

0.4107 

0.25 

0.5210 

0.30 

0.6362 

0.35 

0.7578 

0.40 

0.8865 

0.45 

1.0235 

0.50 

■  1.1702 

0.55 

1.3272 

0.60 

1.4962 

0.65 

1.6780 

0.70 

1.8737 

0.75 

2,0845 

0.80 

2.3120 

0.85 

2.5566 

0.90 

2.8205 

0.95 

3.1039 

1.00 

3.4083 

1.05 

3.7349 

1.10 

4.0852 

1.15 

4.4595 

1.20 

4.8598 

1.25 

5.2869 

1.30 

5.7419 

1.35 

6.2263 

1.40 

6.7408 

1.45 

7.2865 

1.50 

7.8653 

1.55 

8.4773 

1.60 

9.1243 

1.65 

9.8072 

1.70 

10.5270 

1.75 

11.2856 

1.80 

12.0833 

1.85 

12.9214 

1.90 

13.8014 
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0 

.00 

0 

.05 

0 

.10 

0 

.15 

0 

.20 

0 

.25 

0 

.30 

0 

.35 

0 

.i|0 

0 

^5 

0 

50 

0 

55 

0 

60 

0 

65 

0 

70 

0 

75 

0 

80 

0 

85 

0 

90 

0 

95 

1 

00 

1. 

05 

1. 

10 

1. 

15 

1. 

20 

1. 

25 

1.30 
1.35 
1.^0 
1.^5 

1.50 
1.55 
1.60 
1.65 
1.70 
1.75 
1.80 
1.85 
1.90 
1.95 


Table  4.2  -  Continued 

k  =  1 


X 

"k 

\-M^- 

0 

.0004 

0 

.0668 

0 

.1337 

0 

.2006 

0 

.2688 

0 

.3373 

0 

.4067 

0 

.4770 

0 

.5487 

0 

.6213 

0 

.6962 

0 

.7721 

0 

.8497 

0 

.9295 

1 

.0107 

1 

.0942 

1 

.1797 

1 

2669 

1 

.3565 

1 

4484 

1 

5425 

1 

6386 

1 

7374 

1 

8384 

1 

9414 

2 

0472 

2 

1549 

2 

2650 

2. 

3775 

2. 

4926 

2. 

6099 

■  2. 

7299 

2. 

8520 

2. 

9762 

3. 

1031 

3. 

2324 

3. 

3635 

3. 

4976 

3. 

6336 

3. 

7720 

FIRST  DIFFERENCE 


0 

.0004 

0 

.0664 

0 

.0669 

0 

.0669 

0 

.0682 

0 

.0685 

0 

.0694 

0 

.0704 

0 

.0717 

0 

.0728 

0 

.0746 

0 

.0760 

0 

.0776 

0 

.0797 

0 

.0812 

0 

.0835 

0 

0855 

0 

.0872 

0 

.0896 

0 

0919 

0 

0941 

0 

0960 

0 

0989 

0 

1010 

0 

1030 

0 

1058 

0 

1077 

0. 

1101 

0. 

1125 

0 

1152 

0. 

1173 

0. 

1200 

0. 

1221 

0. 

1242 

0. 

1269 

0. 

1293 

0. 

1312 

0. 

1341 

0. 

1359 

0. 

1384 
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Table  4.2  -  Continued 
k  =  2 


X 

\\-^l^^^ 

0.00 

0.0004 

0.05 

0-.0599 

0.10 

0.1204 

0.15 

0.1803 

0.20 

0.2411 

0.25 

0.3019 

0.30 

0.3632 

0.35 

0.4253 

0.40 

0.4877 

0.45 

0.5512 

0.50 

0.6151 

0.55 

0.6799 

0.60 

0.7453 

0.65 

0.8118 

0.70 

0.8790 

0.75 

0.9476 

0.80 

1.0171 

0.85 

1.0875 

0.90 

1.1586 

0.95 

1.2310 

1.00 

1.3043 

1.05 

1.3789 

1.10 

1.4540 

1.15 

1.5308 

1.20 

1.6084 

1.25 

1.6870 

1.30 

1.7667 

1.35 

1.8475 

1.40 

1.9294 

1.45 

2.0122 

1.50 

2.0959 

1.55 

2.1808 

1.60 

2.2666 

1.65 

2.3533 

1.70 

2.4414 

1.75 

2.5302 

1.80 

2.6198 

1.85 

2.7110 

1.90 

2.8025 

1.95 

2.8952 

FIRST  DIFFERENCE 

0.0004 
0.0595 
0.0605 
0.0600 
0.06.08 
0.0608 
0.0613 
0.0621 
0.0724 
0.0634 

0.0639 

0.0648 

0.0654 

0.0666 

0.0672 

0.0685 

0,0695 

0.0704 

0.0712  ■ 

0.0723 

0.0733 
0,0746 
0.0751 
0.0768 
0.0776 
0.0786 
0.0797 
0.0807 
0.0819 
0.0829 

0.0837 
0.0848 
0.0858 
0.0867 
0.0882 
0.0888 
0.0896 
0.0911 
0.0916 
0.0926 
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Table  4.2  -  Continued 
k  =  3 


X 

\\.l^^^ 

FIRST  DIFFERENCE 

0.00 

0.0004 

0,0004 

0.05 

0.0569 

0.0565 

0.10 

0.1143 

0.0573 

0.15 

0.1716 

0.0573 

0.20 

0.2291 

0.0575 

0.25 

0.2873 

0.0582 

0.30 

.  0.3451 

0.0578 

0.35 

0.4037 

0.0587 

0.40 

0.4624 

0.0587 

0.45 

0.5215 

.   0.0592 

0.50 

0.5815 

0.0600 

0.55 

0.6418 

0.0603 

0.60 

0.7026 

0.0608 

0.65 

0.7639 

0.0613 

0.70 

0.8260 

0.0621 

0.75 

0.8886 

0.0626 

0.80 

0.9520 

0.0634 

0.85 

1.0158 

0.0638 

0.90 

1.0806 

0.0648 

0.95 

1.1456 

0.0651 

1.00 

1.2117 

0.0661 

1.05 

1.2781 

0.0664 

1.10 

1.3453 

O.O672 

1.15 

1.4130 

0.0677 

1.20 

1.4815 

0.0685 

1.25 

^ 1.5511 

0.0695 

1.30 

1.6210 

0.0699 

1.35 

1.6913 

0.0704 

1.40 

1.7625 

0.0712 

1.45 

1.8342 

0.0717 

1.50 

1.9064 

0.0722 

1.55 

1.9798 

0.0733 

1.60 

2.0531 

0.0733 

1.65  - 

2.1274 

0.0743 

1.70 

2.2026 

0.0752 

1.75 

2.2778 

0.0752 

1.80 

2.3538 

0.0760            r 

1.85 

2.4308 

0.0770 

1.90 

2.5080 

0.0772 

1.95 

2.5862 

0.0783 
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0 

.00 

0 

.05 

0 

.10 

0 

.15 

0 

.20 

0 

.25 

0 

.30 

0 

.35 

0 

.^0 

0 

^5 

0 

50 

0 

55 

0 

60 

0 

65 

0 

70 

0 

75 

0 

80 

0, 

85 

0 

90 

0 

95 

1. 

00 

1. 

05 

1. 

10 

1. 

15 

1. 

20 

1.25 
1.30 
1.35 
1.^0 
1.^5 

1.50 
1.55 
1.60 
1.65 
1.70 
1.75 
1.80 
1.85 
1.90 
1.95 


Table  ■4.2  -  Continued 


\ 

\.l^-^ 

0 

.0004 

0 

.0556 

0 

.1113 

0 

.1670 

0 

.2225 

0 

.2787 

0 

.3352 

0 

.3917 

0 

.4482 

0 

.5056 

0 

.5629 

0 

.6210 

0 

.6794 

0 

.7380 

0 

.7972 

0 

.8566 

0 

9169 

0 

.9774 

1 

0387 

1 

0998 

1 

1619 

1 

2245 

1 

2875 

1 

3509 

1 

4148 

1 

4791 

1 

5442 

1. 

6097 

1 

6758 

1 

7415 

1. 

8083 

1. 

8755 

1. 

9434 

2. 

0107 

2. 

0786 

2. 

1491 

2. 

2164 

2. 

2867 

2. 

3594 

2. 

4273 

FIRST  DIFFERENCE 


0 

.0004 

0 

.0552 

0 

.0557 

0 

.0557 

0 

.0555 

0 

.0562 

0 

.0565 

0 

.0565 

0 

.0565 

0 

.0573 

0 

.0573 

0 

.0582 

0 

.0583 

0 

.0587 

0 

.0592 

0 

0595 

0 

.0603 

0 

0605 

0 

0613 

0 

0611 

0 

0621 

0 

0626 

0 

0629 

0 

0634 

0 

0689 

0 

0643 

0 

0651 

0. 

0656 

0 

0661 

0. 

0657 

0. 

0668 

0. 

0672 

0. 

0679 

0. 

0673 

0. 

0680 

0. 

0705 

0. 

0673 

0. 

0703 

0. 

0727 

0. 

0680 
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Table  H.2    -    Continued 
k=  5 


X 

\\fl^^) 

0.00 

0.0004 

0.05 

0.0548 

0.10 

0.1092 

0.15 

0.1638 

0.20 

0.2187 

0.25 

0.2734 

0.30 

0.3286 

0.35 

0.3840 

0.40 

0.4395 

0.45 

0.4952 

0.50 

0.5512 

0.55 

0.6077 

0.60 

0.6647 

0,65 

0.7215 

0.70 

0.7791 

0.75 

0.8367 

0.80 

0.8950 

0.85 

0.9537 

0.90 

1.0123 

0.95, 

1.0718 

1.00 

1.1315 

1.05 

1.1914 

1.10 

1.2517 

1.15 

1.3125 

1.20 

1.3738 

1.25 

1.4351 

1.30 

1.4972 

1.35 

1.5596 

1.40 

1.6220 

1.45 

1.6850 

1.50 

1.7486 

1.55 

^   1.8142 

1.60 

1.8761 

1.65 

1.9378 

1.70 

2.0012 

1.75 

2.0738 

1.80 

2.1336 

1.85 

2.2013 

1.90 

2.2692 

1.95 

2.2956 

FIRST  DIFFERENCE 
0.0004 

0.0544 

0.0544 

0.0547 

0.0549 

0.0547 

0.0552  ^ 

0.0554 

0.0555 

0.0557 

0.0560 
0.0565 
0.0570 
0.0568 
0.0575 
0.0576 
0.0583 
0.0587 
0.0587 
0.0595 

0.0597 

0.0600 

0.0602 

0.0608 

0.0613  ■ 

0.0613 

0.0621 

0.0624 
0.0624 
0.0630 

0.0636 
0.0656 
0.0619 
0.0617 
0.0634 
0.0727 
0.0597 
0.0677 
0.0680 
0.0264 
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Since  H^  ^]^_+i   is  convex,  k  =  0,1,2,...  then  H  H   Is 


convex  providing  m>k.   This  follows  by  writing  H  *H 

*=  k  ra 


and  noting  that  the  composition  of  a  convex  Increasing 
function  with  another  convex  function  Is  Itself  convex, 

^ .  5   General  Results  Based  on  Intuition 

An  Intuitive  generalization  of  (4.3.1)  Is  the 

following:   Let  (F  }   „  be  a  sequence  of  cdf's  and 
J   y  £ " 

suppose  that  F    <  P    whenever  y,  <  y^ .   Let  Y-,  and 

^1  (p)  ^2  ^2         1 

Y^   be  r.v.'s  taking  values  In  Q.      Then, 


E     F    ^  ^Y  ^y  ■   (^-5.1; 

^1  n  (p)   ^2  ^2 


providing  Y^f Y^.   That  (4.3.1)  Is  a  special  case  of 
this  Is  clear  upon  letting  Y  =1  w.p.l,  0  < Y,  < 1  w.p.l 
and  F^(x)  =  F^(xy). 

Letting  X^  ^  be  lid  pre-t  r.v.'s  on  v  d.f., 
l  =  l,2,...,n,  (4.5.1)  points  out  how  to  choose  a  set 

of  weights  {a.  }.^,  which  will  maximize  the  pov;er  of 

n 

-=  I  a.X    .   Using  the  notations  of  §4.4,  we  obtain 
1=1     ' 


n 


I  ^'.H   <  I   d  H  providing  I   c.  >    [  d..     (4.5.2' 
'=1  ^  ^  (p)  j=l  -^  '^  j=l  J   i=i  J 


for  s  =  1,2, . . . ,n  and  {c. },{d. }  chosen  to  satisfy  c.  >  0, 

n         n  "^     "^  -^ 

d  >0  and  I    c     =1=     Id.      This  is  seen  by  letting 
'^                   j-1  ^             j=l  J 
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Y^  =  J  w.p.  c.  and  Yp  =  j  w.p.  d   j  = 1,2, . ..,n  ,  in 

^  '^s      s 

(4.5.1).   Clearly,  Pr{Y  <  s}  =  lc.>     I   d .  =  ?v{Y^  <  s] , 

3  J=l  -^   j  =  l  ^       '^ 

so  that  Y^  ^Y^,   As  seen  in  §4.4,  H.   <  H,  whenever 

(P) 
1  < k.   (4.5.1)  now  becomes  (4.5.2). 


It  is  observed  from  Table  3.4  that  if  v  =v  =v, 

s 
then  for  s  = 0,1, 2 , . . . , v  ,  the  quantities  I    n .  are 

j  =  l  ^ 
non-increasing  as  a  approaches  h.      This  would  indicate 

that  the  maximum  power  is  attained  when  equal  weights 

are  assigned  to  each  X    in  T.   That  is,  we  would 

V  ,  1 

choose  Sj  =  ->  j=l,2,.,,,n.   A  preliminary  investiga- 
tion of  the  power  function  of  T  corroborates  this 
finding.   In  a  .05  size  test  the  power  of  T  for  n=2 
is  computed  for  alternative  values  ranging  from  .00 
to  2.00  and  tabulated  in  Table  4.3.   It  is.  seen  that 
when  the  degrees  of  freedom  are  equal,  maximum  power 
occurs  when  the  weights  of  the  pre-t  variables  are 
equal  (a  =  h).      V/hen  d.f.  differ,  the  power  will  be 
maximized  at  some  point  for  which  the  pre-t  variable 
with  the  larger  d.f.  has  the  larger  weight.   It  will 
be  shown  subsequently  that  that  point  depends  not  only 
on  the  d.f.  but  also  on  the  size  of  the  test. 
Notice  that  if  i  <  k  and 

1    if  j  =  i        I  1    if  j  =  k 


^    '0    if  J  ^1     Jo    if  J  ?^k 

in  (4.5.2),  we  obtain  the  result  H.   <  H   for  1  < k, 

^  (P)  ^ 
which  was  shown  in  §4.4. 
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Unfortunately,  (4,5.1)  has  already  been  shown  not 
to  hold  in  general  by  the  counterexample  to  (■4,3.1). 
It  is  of  Interest,  however,  to  Investigate  what  addi- 
tional conditions  need  be  invoked  to  ensure  (4.5.1) 
(and  hence,  (4.3.1)),   This  is  done  by  example  in  §4.7, 


4.6   Prooertles  of  the  Relation    /  v 
' (p) 


Let  F  and  G  be  cdf's  such  that  F  <  G,  as  defined 

(p) 
in  Theorem  4,2,   Thus,  for  any  ae  (0,1),  there  exists 

constants  a  and  3  satisfying 

F(a)  =  G(3)  =  a  , 
and 

1^  >  1  ■     ■' 

93  " 

We  look  at  some  properties  of  the  relation  "  <  ": 

(P) 
For  cdf's  F,  G  and  H,  let  P  <  G  and  G  <  H. 

(p)        (P) 
Then  F  <  H  (transitivity  of   <  ) ,   This  is  obvious 

(p)  (p) 

by  letting  a,  3  and  y   be  such  that  F(a)  =G(3)  =H(y)5 

so  that 

3a  _  5a  93  ^  . 
9y    93*9y  " 

Since 

F(a)  =  G(3)  =>      F(a-d)  >  G(3-d),Vd>0  (4,6.1) 

is  equivalent  to 


F  <  G  ,  (4.6.2) 

(P) 
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v<je  may  write  for  any  positive  constant  c, 

cF  <  cG  (4.6.3) 

(P) 

to  mean 

cF(a)  =  cG(3)  =^cF(a-d)>cG(3-d),   Vd>0 

which  is  equivalent  to  (4.6.1).   Thus  ( ^1 .  6  .  2 )  and  (4.6.3) 

have  Identical  meanings.   Assuming  also  that  we  may  add 

the  same  Quantity  to  both  sides  of   <  ,  we  have 

(P) 

F<G=^F  +  H<G+H  (4.6.4) 

(p)  (P) 

and  the  following  properties  hold: 


(i)   F.G.,  1  =  1,2  ^  F^+F   <  G^+G^. 
(P)  (p) 

proof: 

F^  <  G^  =^  F^+F^  <  Gi  +^2     ^^    (4.6.iJ) 
(P)  (p) 

F^  <  G^  =^  G  +F2  <  G  H-Gp     by  (4.6.4) 
(p)  (p) 

•'•  F,  +  F„  <  G,  +  G^        by  transitivity. 
^    "^  (P)  ^    ^ 


(11)   For  0  <  d  <  c  <  1, 

F  <  G  ==>cF+(l-c)G  <  dF+(l-d)G 
(P)  (P) 

proof: 

(c  -  d)  >  0   => 

(c-d)F  <  (c-d)G  =  (l-d)G-(l-c)G 
(P) 
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•••  cF+(l-c)G  <  dF+(l-d)G   after  2  apolicatlons 
(p)  of  (^.6.4) 

(ill)   Let  {c^},{d^},  i  =  l,2,.,.,n  satisfy 

n      n 
c>0,  d.  >0,   Ic.=  I   6.  =  1,    and  for  some  k, 
1=1  ^   1=1  ^ 

c^  >  d.    for  1  =  1,2,  . . . ,k 
while 

c.  <  d.    for  1 =  k+1, . . . ,n. ' 

If  F.   <  F.  for  1  <j  then  • 


11  11 

y  c.F.   <    y  d.F. 

.=1     (p)  1=1 


proof : 


For  1  =  1^2,. ..,k,  (c.-d.)>Oso  that 


(c,-d.)F.  ^.^(c.-d.)Fj^^^.        (il.6.5) 


For  1  =  k+2,  ...,n,  (d.  -c.)>Oso  that 


(d,-c,)F^^^   <  (d.-c.)F.  .        (11.6.6) 
(P) 

Applying  property  (1)  to  (11.6.5)  k-1  times  and 

to  (4.6.6)  n-k-2  times,  we  obtain  . 

k  k 

I    (c     -d    )F       <       I    (c.-d    )F  ,       (4.6.7) 

1=1      ^        ^      "-  (p)  1  =  1      1        1      i^+J- 

and 

I      (d^-c    )F  <         I      (d    -c    )F        (4.6.8) 

l  =  k+2      ^         ^      ^^^  (p)  l=k+2      111 
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Applying  property  (1)  again  to  (i|.6.7)  and  (^.6.8), 

I    (c.  -d  )F   +    I      (d  -c  )F 
i  =  l   ^    11    1=1^+2   1    1   ^+i 


<    I  (c.  -d  )P     +    I      (d  -c  )F 
p)  i=l   ^    ^   ^^    irk+2  ^    ^   ^ 


which,  by  (^.6.^),  becomes 

k  n  n 

I    c.F.  +       I      c.F.  -       y      c.F,  ^, 
i=l    ^    ^      i=k+2    ^    ^        1^1      1    ^+1 

i7^k+l 

(^.6.9) 
k  n  n 

<         I    d.F.  +       I      d.F.  -       y      d.F,  ^,. 
(p)    i  =  l    ^    ^      i=k+2    ^    ^        1  =  1      ^    ^-^1 

l7^k+l 


Noting  that 


^k+1   ^   1  - 

1 


.^,^1'      \+l   =    ^-     .^/l' 
1=1  1=1 

i?^k+l  i7^k+l 


(^.6.9)  gives  the  result. 

^.7   Some  Counterintuitive  Examples 

From  Table  3.^  it  is  seen  that  if  v^  =Vp  =  v,  the 
weights  ^^'^j_(o')  ^  j^  =  o  corresponding  to  the  linear  combina- 
tion T  =  aX   +(l-a)X   are  such  that,  for  some  k, 
^1         ^2 

n^(a)  >  n^ih)  ,       i  =  0,1, . . . ,k  , 
and 

ri^(a)<n^('^),   i  =  k+l,k+2,  .  .  .  ,v  . 
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Thus  by  the  results  of  §^4.4  and  property  (111)  of  §^.6, 
we  would  have 


V  V 

I   n.(a)H   <   [  n.(is)H   ,         (4.7.1) 
1  =  1  ^    ^  (p)  1  =  1  ^    ^ 


which  is  to  say  maximum  power  occurs  v;hen  a  =  h.      An 

n 
extension  to  T  =  I   ^ •  ^v)   would  suggest  a  choice  of 

1  1=1  ^   1 

a .  =  — ,  1  =  1 , 2 , . . . , n  whenever  v .  =  v ,  V 1 . 

These  results  are  very  concordant  with  our  Intui- 
tion.  It  seems  reasonable  to  weight  equally  those 
variables  carrying  equal  amounts  of  information  and  to 
weight  more  heavily  those  which  carry  a  greater  store 
of  Information.   The  result  ( iJ .  7  . 1 )  rests  on  the  assump- 
tion {H.6.H),    which  itself  appears  very  much  in  agree- 
ment with  intuition.   If  sampling  from  the  distribution 
G  results  in  a  more  powerful  test  than  does  sampling 
from  the  distribution  P,  the  natural  inclination  would 
be  to  sample  from  a  mixture  pG  +  (1  -p)H  given  a  choice 
between  it  and  the  mixture  pF  +  (l-p)H.   Unfortunately, 
this  is  not  always  the  best  thing  to  do;  i.e.,  (4.6.^4) 
is  not  true  in  general.   Although  a  computer  study 
Indicates  (^1.7.1)  is  still  valid,  the  extra  conditions 
which  guarantee  the  validity  of  .(4.6.4)  are  elusive. 
The  remainder  of  this  section  is  an  exploration  of 
some  counterexamples  to  (4.6.4). 
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Counterexamples  to  (^  .S.k): 


Let  F  and  G  be  cdf's  and  F  <  G.   The  assumption 

(p) 
of  (^.6. -4)  Is  for  any  cdf  H,  we  may  declare  P  +  H  < 

(P) 
G  +  H.   The  examples  which  follow  take  K  as  either 

F  or  G. 

(1)      Let    G   be   the   cdf   of  a  r.v.    whose   density   is 

.  ■^c"     ,  I  X  [  <  c 

:(x)    =   -!  .Ic"^  ,       e  <  ]x|  <  2c 
0        ,        otherwise 

and   let   F(x)=G(|).      Choose   x   and   y   such   that 
F(x)=G(y)=a.       For    d>0,    F(x-d)-  G(^^-=-^)    > 

G(|- d)  =  G(y  -  d),    so   that   P     <     G.      Let   a,    a   be 

(P) 
such    that 

^/2{F(a)  +G(a)}    =   is{G(3)  +G(3)}    =    .95. 
We   have 

G(e)  =  .95     =^    ^  =  |c 
i§{F(a)  +  G(a)}  =  .95     =^    a  =  2c. 

But    for   d  =  .2,    c  >  .2,    we    have 

%{F(a  -d)  +G(a  -d)}    =   Jg{G(c  -  .  1 )+  G  (2c  -  .  2  )  } 
=    .95  -— ^<  .95  -^=G(1.5c  -  .2)  =G(3-d) 
so   it    is   not    true   that 


P<G=*>      F  +  G<G+G 
(P)  (P) 
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(ii)   Let  G  be  the  cdf  of  the  uniform  r.v,  with  density 

r 

S  ,    -1  <  X  <  1 

:(x)  = 

0  ,   otherwise 


and  let  $  be  the  cdf  of  the  standard  normal  r.v, 
Choose  X  and  y  such  that  $(x)  =G(y)  =o.      Letting 


dx 


dx  ^  g(y)  ^  /2lf  e  ^^ 
dy    (t)(x)  ■     2 


>  1 


so  that 

0  <  G  ,  ■ 
(p) 

Let  a  and  3  be  such  that 

^^{<I>(a)  +<D(a)}  =  ^-2{f  (3)  +G(3)}  =  .95 


If  Z   represents  a  number  such  that  $(Z  )=l-p, 
we  have 


and. 


a  =  Z_05  ^  ^-^^5 


=  Z^^Q  =  1.281. 


For  several  values  of  d,  ^(a-d)  and  ^^{^(g-d) 
+  G(3-d)}  were  calculated  as  follovjs  : 


05 

.9446 

10 

.9388 

15 

.9326 

20 

.9258 

281 

.9137 

(g-d)     ^-2{J'(g-d)+G(3-d)} 

.9454 
.9406 
.9352  . 
.9300 
.9206 


It  is  seen  that 

$(a-d)<?5($(3-d)+G(6-d)} 

for  the  cases  considered,  so  it  is  not  true  that 

$  <  G  =>$+<|)  <  G  +  f 
(P)  (p) 

Interestingly  enough,  however.  It  is_  true  that 

$  <  G  =^$+G  <  G  +  G. 
(P)  (p) 

To  see  this,  let  a,  g  be  such  that 

Jg{$(a)  +G(a)}    =   G(B)  ■ 
Then, 


3a  _    2g(g: 


>  1 


3l  "  4)(a)  +  g(a) 

(ill)   Let  G  be  any  cdf  satisfying: 

G(-c)=0,  G(c)=l  and  F  <  G  , 

(P) 

where  F  is  the  cdf  of  a  normal  r.v.  with  mean  0 
and  variance  o^ .      Let  a   and  3  be  such  that 

Jg{F(a)  +P(a)}  =  i5{F(3)  +G(B)}  =  1  -  *  .  , 
where  (J)  is  chosen  so  that  6>c. 

F(a)  =  1-4)  =t>   a  =  aZ^ 

^-slFCB)  +  G(3)}  =  1  -  (1)  =^   3  =  aZ^^ 

2<p 


Now,  letting  f(x)  =  ^^^  , 

ff^^.l   Iff   f(3).2f(a) 


iff   e  ^^"  >  2e 


a 
2o' 


iff   a^  -  6^  >  20^  In  2 


iff   Z^  -  Z^  ,  s  2  In  2  =  1.3863. 


2    2 
For  several  values  of  <^ ,    Z  -  Z„ ,  was  calculated 


as  follows 


Z 


'?f 


2    2 

_^> ^ 


.20      .8i]2      .253  .645 

.10     1.282      .842  .935 

.05     1.645     1.282  1.063 

.025    1.960     1.645  1.136 

.01          2.327          2.054  1.196 

.005    2.575     2.327  1.216 

Thus  for  the  cases  indicated,  ~<1,    so  it  is 

0  P 

not  true  that 


P  <  G 
(P) 


F  +  F  <  F  +  G 
(P) 


(iv)   Let  $  be  the  cdf  of  the  standard  normal  r.v. 

1  ? 

with  density  (p(x)  = exp(-hx    )   ,  and  G  the  cdf 

/2?  2 

of  a  normal  r.v.  with  density  g(x)  = expC-T^-^-), 

/2¥y      ^^ 


$(x)  =  G(yx) 


For  0  <  Y  <  1,  $  ^  G  since  if  x,  y  are  such  that 
$(x)=G(y),  then  y-YX  or  ^  =  ->l. 
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Let  a  and  3  be  such  that 


J2{<J>(a)  +  $(a)}  =  hi^iB)  +G(B)}    =   1- 


Then,  a  =  Z^  and  since  n(()(nx)^0  as  n->o°,  y  may 
be  chosen  small  enough  so  that  -  (j)  (-)=  0 .  00000 
and  <I'(:^)  =  G(3)  =  1.00000,  we  have 


-  ^2      (to  5  decimal  places 


3B       2(t)(a)        2(^(a)        2({,\a) 


For  several  values  of  p,  ||-  was  computed  as 
follows: 


p 

*(Zp) 

*.(z,p) 

da 
33 

.200 

.2799 

.3864 

69 

.100 

.1754 

.2799 

80 

.050 

.1031 

.1754 

85 

.025 

.0584 

.1031 

88 

.010 

.0266 

.0484 

91 

.005 

.0145 

.0266 

92 

For  each  of  these  cases  we  see  |^  <  1 ,  so  it  is 
not  true  that 


<  G 

(P) 


<  $  +  G 
(P) 


Like  example  (11),  however,  we  note  that  it  is 
true  that  for  appropriately  chosen  y , 


<     G 
(P) 


$  +  G   <  G  +  G 
(P) 


To  see  this,  let  a  and  3  be  such  that 


hi^ia)   +  G(a)}  =  G(3)  =  1  - 


so  that  3 =  yZp  and,  if  y  is  chosen  so  that 
~<^  {-)  =   0.00000  and  $(^)  =  1.00000, 


a  =  Z 


2p 


In  this  case 
8a 


2M)  2    ^'^^^p)  ^  . 

(a)  +g(a)    Y  ^(Z^p)  -   ' 


Vp 


(v)   Let  H^  be  the  cdf  of  X2^^^  for  1  e  {j,k}.   If 


j  <  k,  we  have 

H.   <  H 


k 


by  §4.4 


^   (p) 

Let  j=2,  k = 48  and  a  and  3  such  that 


^2{H  (a)  +H  (a)}  =  ^-^{H  (B)  +  H,  (3)}  =  .90 


We  find 


and 


a  =  .65994 
6  =  .41142 


3^-^2(g^-^^8^g^. 53092  +  .00185    Qpo,,  ' 
83   "   2h2(a)     "   2(. 28694) .9284<1 


So  it  is  not  true  in  general  that 


H.  <  H, 
'   (P)  ^ 


H.  +H,   <  H.  +  H, 
J    J  (p)  J    k 


8^ 


^.8   The  Power  of  aX   +(l-a)X 

^1         ^2 


Write 


T   =  ex^,  +  (1  -  9)X      for   e  e  (0,1)  . 
Using  §2.1,  we  get 
EiTl"")    =  (4.8.1) 


If  v^y^v^    this  quantity  is  minimized  by  a  value  of  6, 
say  6^,  which  depends  on  n.   For  example, 

^1-2 


1    v^  +  v^  -  4  ' 

while  Q^,d    ,,..    are  more  complicated  functions  of  v 

and  v^'      Thus,  there  exists  no  unique  value  of  0  which 

2n 
minimises  E(Tg  )  for  all  values  of  n  (1,2,3,...).   From 

Theorem  4.3,  then,  it  follows  that  there  is  no  linear 

combination  of  pre-t  r.v.'s,  Tg,,  which  has  greater  power 

than  all  other  linear  combinations  for  all  size  tests. 

If  ^2=^2  '^his  ^^y   1^0  longer  be  the  case.   The 

value  of  0  which  minimizes  (4.8.1)  is  found  by  solving 


^EiTl^) 


=    0 


Letting  V  be  the  common  value  of  v   and  Vp ,  we  have 


Calling  the  product  of  the  first  five  factors  in  the 
summand  G^  and  letting  S^  =  (i-^),  we  notice  that 


and 


G.  =  G   . 
1    n-i 


S.  =  -S   . 
1     n-i 


so  that  9=3^  is  a  root.   There  can  exist  no  other  roots 
by  the  symmetry  of  G.. 

'?r\ 

Thus  for  all  n  (=1,2,3,...).  E (T^    )    is  minimized 
by  e^h.      Thus,  the  test  Tj  is  admissible  when  v,  =  v„ 
since,  by  Theorem  ^.3,  there  can  exist  no  test  T„ , 
(6'  T^Jg)  more  powerful  for  all  sizes. 


CHAPTER  5 
COMPARISONS  WITH  OTHER  METHODS 


5-1   Ghosh's  Results 

Ghosh  [7]  expressed  the  cdf  of  the  sum  of  two 
lid  t-variables,  say  t^   and  t^,  as  a  finite  sum  of 
hypergeometrlc  functions  v:hlch  reduces  to  explicit 
formulae  for  odd  degrees  of  freedom.   To  do  this  he 
v;rltes  the  density  of  t'  +t"  as 

V     V 


f(t) 


f^(x)f^(t-x)dx 


where  f^,(y)  Is  the  density  of  t  .  After  a  formal 
series  expansion,  a  term  by  term  Integration,  and 
use  of  the  duplication  formula,  this  becomes 

is(2v  +  l) 


f(t)  = 


^v 


^{^}^{2^^} 


Hfi  2v+l  v+2 


Hv+t' 


(5.1.1)' 


where,   H(a,b,c;  x)  = 


r{c} 


I 


r{a+i}r{b+i}  1 


rTaTNbJ^^Q     riTTc+iT 


Is  the  hypergeometrlc  function.   Integrating  (5.1.1) 
the  cdf  Is  expressible  as 
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Pr{t'  +  t  <  t} 

V     V 


1    r{^} 

2'rv — "— 


r(2i±l}r{2::^±|i±i} 


2^/?r2{^}i^o   i!r{^^i2i+2} 


t2 


p- 


(5.1. 


l-is 


(l+^)M2v  +  21  +  l) 


■du 


for  t  >  0.   If  V  is  odd  this  reduces  to  an  expression 


1  1  3 


111 


involving  U{^,j,^;    x)  or  H(J,J,|;  x)  wher 


t' 


4v+f 


If  V  is  even  the  hypergeometric  functions  reduce  to 

expressions  involving  the  complete  elliptic  integrals 

of  the  first  and  second  kinds  given  by 

Tr/2 


K(x) 


E(x)    = 


2 

(1- X  sin    I 

0 

7T/2 

p  '2 

(1- X  sin    6)    d3 


d3        and 


In  particular. 


Pr{t'  +t''  <  t}    =   k  +  -  tan~^(S 


■1        1 


2  TT 


(^) 


Pr{t2  +tl<t]    =    ^  +   ^{(i+q)E(q)  -  (l-q)K(q)} 

p^r^'    ,   ^-"   ...     _    1     ,     2/3t(l8+t^)    ^    1    .        -1/    t    , 
Frit,  +  t^  <  t)    =   ^  +   — =i-^ M      +   -  tan      (0-7^5:) 


I       .  'I 


Pr{t.  +t^-  <  t)    =   ^  +  — -L^  {(8p"-31p^+48p^  +  5p+2)E(p) 


'or    t  >  0   and   q=^-p    ,    P=— ^ — ^ 
8+t  16+t 


-(^p   -l6p3+6p'^  +  i|p  +  2)K(p)} 
2 


Ghosh  gave  tables  of  the  cdf  of  t'  +t"  using 
these  formulae.   For  v>^,    he  uses  numerical  Integra- 
tion to  evaluate  (5.1.2).   His  values  agree  to  four 
decimal  places  with  those  found  using  the  method  of 
Chapter  3.   To  approximate  the  distribution  of  t'  +t" 
v;hen  v^  <  v^ ,  he  expands  the  density  f(t)  to  order  v""^ 
and  integrates  term  by  term.   For  the  few  cases  he 
tabulated  there  appears  to  be  fairly  good  agreement 
with  the  exact  values. 

5 . 2   The  Normal  Approximation 

Because  t^  has  normal  limiting  distribution,  so 
does  the  statistic  T  given  by  (3-1.2).   Of  interest  is 
how  large  the  Vj ' s  must  be  to  adequately  approximate 
the  distribution  of  T  using  the  normal  distribution. 
For  n=2,  (3.1.2)  becomes 

T  =  aX    +  (l-a)X^,  (5-2.1) 

1         ^2  .  ■       ' 

and  the  quantity 

V  =  -^  +  -^^^  (5.2.2) 

represents  the  variance  of  T.   Letting  Y  be  a  normal 
r.v.  v;ith  zero  mean  and  variance  v,  the  cdf's  of  T 
and  Y  are  compared  for  several  values  of  a,v,  and  Vp. 
It  was  found  that  the  best  approximation  occurred  when 
v^   =   v^    =   \)   and  a  =  ^  -   In  this  case  (5.2.2)  becomes 
(2v-^)"  ,  and  the  com.parison  is  given  in  Table  5.1. 
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The  maximum  error  was  .066  for  v=3,  .017  for  v=5 
and- .010  for  v=7,9  or  11.   When  v=13  or  15,  the 
largest  absolute  deviation  is  seen  to  be  .006.   The 
approximation  becomes  somewhat  better  in  the  tail  areas 
with  an  average  error  of  less  than  .002  forv>9.   The 
normal  approximation  appears  adequate  to  two  decimal 
places  V7hen  v^l5.   This  may  not  be  the  case,  however, - 
if  the  weights  or  d.f.  .differ. 
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Table  5.1 
COMPARISON  OF  THE  CDF  OF  igX'  +  i-iX"  V/ITH  A  NORMAL  CDF 

P^(x)  =  cdf  of  isX^  +isX^ 
P2(x)=cdf  of  Y  ~ N(0,{2v-4}~  ) 


.1 

3 

.579 

.556 

11 

.671 

664 

.2 

.654 

.611 

.810 

802 

.3  ■ 

.721 

.664 

.903 

898 

.4 

.778 

.714 

.956 

955 

.5 

.826 

.760 

.982 

983 

.6 

.864 

.802 

.993 

995 

.7 

.894 

.839 

.997 

999 

.8 

.917 

.871 

.999     1 

000 

.9 

.935 

.899 

1.000 

.1 

5 

.606 

.597 

13 

.686 

681 

.2 

.704 

.688 

.832 

826 

•  3 

.786 

.769 

.923 

920 

.ij 

.850 

836 

.969 

970 

.5 

898 

890 

.989 

991 

.6 

931 

929 

.996 

998 

•  7 

954 

957 

.999   1. 

000 

.8 

969 

975 

1.000 

•  9 

979 

986 

.1 

7 

629 

624 

15 

.700 

695 

.2 

743 

737 

.851 

846 

.3 

833 

829 

.938 

9'37 

.4 

897 

897 

.978 

979 

.5 

939 

943 

.993 

995 

.6 

965 

971 

.998 

999 

.7 

980 

987 

1.000    1. 

000 

.8 

988 

994 

.9 

993 

998 

.1 

9 

654 

646 

.2 

783 

773 

.3 

877 

869 

.4 

936 

933 

.5 

969 

969 

.6 

985 

988 

.7 

993 

996 

.8 

997 

999 

.9 

999   1. 

000 

Chapter  6 
APPLICATIONS 

6 . 1   The  Behrens-Fisher  Problem 

Using  the  sampling  scheme  for  the  Behrens-Fisher 
problem  described  in  section  1.2,  the  test  statistic  . 
for 

Hq:  d  =  0 
vs  .     H-,  :  d  >  0  , 

where  d  is  the  difference  in  normal  population  means. 


y-|_-y25  IS 


N. 


;i 


"  =  ^M,"ij^ij-J/ij'ij^  •         ^'-'-'^ 


To  find  the  distribution  of  T,  we  use  the  property 
that  if  R  -  Xy  and,  conditional  on  R,  U~N(0,R"-^),  then 
v^U~t^.   This  may  be  demonstrated  as  follows: 


Pr{A7U  <  x}    - 


Pr{U  <  —    I    R=r}dPr{R<r} 


v^ 


0 


f/AJ 


2 
pU    r 


^-sv-l 
r  e 


0     -» 


(2-n)'         r{J}2- 


du  dr 
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r  //^   ,-Js(v+i) 


J        r{h}T{^] 


;^(v-i)    -i^r(l+u^)   .     , 
r  e  ^  dr  du 


/a;   r{^^} 


r{J}r{J2} 


'--  (i  +  u2)-^'^^-+i) 


du 


which   is   the   cdf   of  t    //v.      Thus    AJ  U  ~  t      (or  U~X    ) 

V  V  V 

From  §1.2,  we  have 


"^  ~  ^n.-l 
a. 


(6.1.2: 


independently  for  i=l,2,  and,  conditional  on  SS . , 


2 
a . 


^IJi^ij^ij-^i  1~  ^(°'s#r^  ^6-i-3) 


independently  for  i=l,2,  so  that 

.  N. 


-^   I   y  a.  .X.  .  -  u. 


(6.1.i]) 


Independently  for  1  =  1,2,  where  v.  =n.  -1.   Therefore, 
we  may  think  of  T  as  having  the  representation 


T  =  a,X   -  a^X   +- 

1  V.     2  Vp   T 


(6.1.5) 


where  X^  's  have  independent  pre-t  distributions  on 
i 

V.  d.f. 

1 
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To  test 

vs .    H  :  d>0  , 


we  use  T  as  test  criterion,  noting  that  when  H^  is 
true,  T  has  the  same  distribution  as 


a.,  A   +  a^ X 
1  V-,    2  v^ 


so  that  the  methods  of  Chapter  3  are  applicable. 
Notice  that  the  power  function  of  T  is  given  by 


Pr{T>A,}  =  PrCa^X   +  a^X   >A^--} 


2 
which  is  independent  of  the  unknown  variances,  a,  and 

a^ .   The  critical  values,  A  ,  satisfying 


Prla^X   +  a„X  >  x    }    =    ^ 
1  V.    2  ^p    ^ 


are  given  in  Table  3-5  for  4)  =  .10,  .05,  .025,  and  .01 

Yet  to  be  discussed  is  the  selection  of  initial 
sample  sizes  n,  and  n„,  and  the  constants  t  and  {a..} 
j=l,...,N..   If 

SS. 
1 

<  n, 


2  2  ^  '^i 
a .  T 

1 


we  take  N.  =n.  +1  and  if 


we  have 


-j\  >   n 
a.T 
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N.  = 
1 


"SS. 
a_ 

2  2 

a .  T 
-  1  — ' 


+  1 


Let 


2  2 
a .  T 
1 


and  denote  by  xj  a  chl-square  r.v.  with  v  d.f.   Noting 
that 


SS. 

1        2 

2  '^  ^n.-l  ' 
a.      'i 
1 


we  have 


SS. 

V2  <   -^i    iff    X^.-l<  "i^i 
a.T  1 


For  k  >  n.  +1 
1 


SS. 


Pr{N,  =  k)  =  Pr{k  -  1  <  -^  <  k} 


2  2 
a.  T 

1 


=  Pr{(k  -  1)(|)^  <  X^  <k(J)^},   v^=n^-l 


Thus, 


ffCN.  )  =  (n^  +  l)Pr{x^  <  n^<p^} 


+    I 


k==n.+l 

1 


(k-l) 


(6.1.6) 
+  1  1  dF(x') 
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where  F(x^)  is  the  cdf  of  x^  •   The  summation  reduces  to 

i 


ri .  4  . 
1^1 


'SS. 
i_ 

2  2 

a.  T 

-  1  — 


+  1   dF(x') 


which  may  be  represented  as 


x'dP(x')  +  R 


where  R  is  such  that 

0  <  R  <  Pr{n.({).  <  x^  }  <  1 

1 

Neglecting  this  fractional  remainder. 

E(M.)  =    (n.  +l)Pr{x^  <"i*i^  +  F" 

i  i 


xMf(xM.  (6.1.7) 


n .  d) . 
1^1 


The  expected  total  sample  size  N =  N,  + Np  is  smallest 
when  we  minimize  ^(N, )  and  ^(Np).   Since  the  expression 
£■(1^)  depends  on  the  unknown  parameters  o.,  it  is  not 
clear  how  the  initial  sample  sizes,  n.,  should  be  chosen 

to  minimize  ffCN)  when  no  prior  information  is  given 

2 
about  the  o  's.   Seelbinder  [1^4  ]  investigated  a  related 

selection  procedure  and  produced  a  set  of  tables  indi- 
cating the  effect  of  Initial  sample  size  choice  on  the 

2 
total  sample  size  for  various  levels  of  o  .   The  tables 
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are  of  value,  however,  only  when  we  know  the  approximate 

2 
magnitude  of  a  . 

The  constant  t  may  be  appropriately  chosen  to  give 

a  confidence  interval  for  d  =  p  -  y   a  preassigned  width 

say  2v7.   This  is  accomplished  by  selecting  x  to  satisfy 

Pr{-X   <T-f.A   )  =  !-♦, 
2  2 


or. 


N,         Nj 


"''"'*/'  %-^i"ij'iJ  -  J^^aj'^jj  -d  =  %/'  =  !- 


This  statement  reduces  to 


Nj         Ng 


'"-^ jU^iihr /^^W^2i  -  \^j  <-  -3 


SO  that  a  100(1  -cj))^  confidence  interval  for  d  of 
v/ldth  2^fJ   may  be  obtained  by  selecting 


^.  (6.1.9) 

^2 


Finally,  it  should  be  noted  that  the  selection  of 
the  constants  {a^  },  j  =  1 ,  .  .  .  ,  N^^ ,  (i  =  l,2)  to  satisfy 
the  constraints  (1.2.2)  may  be  arbitrary.   To  eliminate 
this  arbitrariness,  one  might  use  the  additional 
condition 
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^^n.+l  =  «i,n.+2  =    '"    ^  '^•i,N.  '       (6.1.10) 
'  1       '  1  '  1 


which  effectively  weights  all  information  in  the 
second-stage  sample  equally.   Now,  if 


a  =  a .  .    for  j  =  1 .  2 ,  .  .  .  ,  n . 
ij  '  '      1 


ana 


b  =^  a  .  .    for  j  =  n  .  +1 ,  n  .  +2  ,  .  .  .  -  N . 


we  have. 


N. 
1 

y  a.  .  =  n.a  +  (N.  -  n. )b  =  1 

j'l  -J  ^  ^        ^ 


N.  2  2 

la..    =n.a+(N.-n.)b   = '  ^^—  =  f~      , 

so  that 

n^  ±  (n^  -n.N^d  -  (N.  -n.)f~^)}''^ 

1  1 

and,  (6.1.11) 

1  -  n.  a 

b  =  r. i-. 

N .  -  n  . 
11 


From  these  solutions  it  is  recommended  to  choose  the 
one  minimizing  |a-b|. 


Numerical  Example: 

An  initial  sampling  of  6  observations  from  each  of 
the  normal  populations  gave  the  results: 

I        ■  II  . 


20.6 

19.7 

19.3 

20.1 

20.2 

19.8 

21.4 

18.9 

20.9 

20.3 

21.7 

18.5 

It  is  desired  to  test  for  a  difference  in  popu- 
lation means  and,  if  significant,  to  estimate  the 
difference  to  within  1  unit  with  probability  .95. 

From  the  data  we  calculate 

6  6 

I   X-L-  =124.1        I   X   =117.3 
j=l  ^  j=l  '^^ 

SS^  =  3.7483        SS2  =  2.4750 
n^  =  6        .       n^  =  6 

In  accordance  with  the  results  of  Chapter  4,  we  will 
choose  a-^^=a2=^.   Computing  t  by  (6.1.9)  ,  since  we 
will  use  this  value  to  obtain  the  confidence  interval 
if  there  is  a  true  difference. 


\025    -^^^ 


1.2330 


where  X  ^^^    is  obtained  from  Table  3.5.   V/e  calculate 
f^  =  9.86216   and   f2  =  6.51192 


where 


SS. 
f.  ^    ^ 


1     2  2 
a.T 


This  gives  N,  =10,  N^  =  7 ,  so  that  four  additional 
observations  are  required  from  the  first  population,, 
while  one  more  is  needed  from  the  second.   Using 
(6.1.11),  we  obtain 


a,  =a„  =  •••  =ag=  .IO965I 


"I7^"l8  =  °'l9  =  "l,10=  -0855211 


and 


Op-,  =  Kpp  =  •••  =cip/-  =  .15882 


a^„  =   .04708  . 

The  results  of  the  second-stage  of  the  sampling  process 

are  : 

I        II     , 

20.7      18.4 

20.3 

19.1 

21.4  ■     . 


To  test 


Hq:  v^  =  V^ 


vs.   H^:  y^>y2 


100 


we  use  as  test  statistic 


.  ,  10         7 


rejecting  H„  at  a  significance  level  (t>    if  T>  X,. 
Computing  T=  .8776  and  using  Table  3-5,  vje  may  reject 
H  at  (j)  =  .025.   From  (6,1.8)  a  95?  confidence  Interval 
for  li-,  -  vip  with  desired  bound  of  1  unit  is 

1.082  ±  1. 

6 . 2   A  Nonstationary  First  Order  Autoregressive  Process 

The  model  for  the  nonstationary  first  order,  auto- 
regressive process  may  be  written  as 


^t  =  \^t-l  +  ^  '      t  =  l,2,. 


^0  ~  ^0 


.  ,m 


(6.2.1) 


Suppose  that 


X  = 


has   the   m.ultivariate   normal    distribution   with    zero   mean 
and   dispersion  V,    and 


e^    ~   N(0  ,  o^ 


independently  for  each  t  and  also  independent  of  X 
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The  subscript  on  a   indicates  the  nonstationarlty  of 

the  process.   One  test  of  interest  is  whether  the  process 

is  indeed  nonstationary .   The  null  hypothesis  would  be 


for  some  constant  n^,  and  the  test  statistic  and  its 
distribution  are  developed  in  [5,12].   Another  hypothesis 
of  interest  is  that  a,  remains  the  same  for  two  partic- 
ular time  periods,  say  the  r    and  the  s   .   The  null 
hypothesis  now  is 


H„:a=a,      l<r<s<m. 
0    r    s 

If  s  =  r  +  1  for  example,  we  are  hypothesizing  that  the 
X  's  have  a  constant  autoregressive  relationship  over 
two  consecutive  time  periods.   The  alternative  hypothesis 
might  be  that  this  effect  is  decreasing: 


H-,  :  a-  >  a 
Irs 


The  experiment  consists  of  taking  replicated 
observations  X,  ^Xpj.-.jX   where  n>m  +  l.   The  repli- 
cations are  necessary  to  make  inferences  regarding 
the  parameters  (ct,  },  k  =  l,...,m.   V/e  calculate  the 


matrix 


n 
W  =  I    XX'  , 


which  may  be  uniquely  decomposed  as 


W  =  GDG' 
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G   vjlll  be  lower-triangular  with  unit  diagonal  elements 

and  D,  diagonal.   These  matrices  will  provide  estimates 

of  the  parameters  a-,,a„j...ja  and  0^,0-,  ....  ,o    .      In 
12,,   'm      O'l'    'm 

particular,  if  g   is  the  i-j    element  of  G,  g   ,  _, 
is  an  unbiased  estimator  of  a^  .   Because  V/  has  the 
Wlshart  density  with  dispersion  V  and  n-1  d.f., 

o 

k     ' 

for  k=l,2,...,m,  independently  (see  [8]);  i.e., 

a. 


independently.   Letting  d  =  a  -  a  ,  we  m.ay  test 


H^:  d=0 


vs.    H^:  d>0 
by 

^  =  ^^^Sr.r-l-gs,s-l)-^^'    -<^ 


where 


-,     a     a 


Under  H„, 


r-1    s-1 


^  ~  ^l^v/^2^V2 
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v/here 


i  €  {r, s}  , 


and  V  =n-i.   V/e  reject  H^  in  a  size  *  test  if  T>X,, 

-1-  U  (p 

using  Table  3-5' 
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